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Abstract 



We study the real-time nonequilibrium dynamics in hot QED plasmas imple- 
menting a dynamical renormalization group and using the hard thermal loop 
(HTL) approximation. The focus is on the study of the relaxation of gauge and 
fermionic mean fields and on the quantum kinetics of the photon and fermion 
distribution functions. For semihard photons of momentum eT <C k <C T we 
find to leading order in the HTL that the gauge mean field relaxes in time 
with a power law as a result of infrared enhancement of the spectral density 
near the Landau damping threshold. The distribution function of semihard 
photons in linear response also relaxes with a power law, with a power that is 
twice that for the mean field. The dynamical renormalization group reveals 
the emergence of detailed balance for microscopic time scales larger than 1/k 
while the rates are still varying with time. The quantum kinetic equation for 
the photon distribution function allows us to study photon production from 
a thermalized quark-gluon plasma (QGP) by off-shell effects. We find that 
for a QGP of temperature T ~ 200 MeV and lifetime 10 ^ t ^$ 50 fm/c the 
hard (k ~ T) photon production from off-shell bremsstrahlung {q — > q'j and 
q — ► qj) at 0(a) grows logarithmically in time and is comparable to that 
produced from on-shell Compton scattering and pair annihilation at 0(aa s ). 
Hard fermion mean fields relax as e~ aTt ln ( w -p*) with up = eT/3 the plasma 
frequency, as a consequence of the emission and absorption of soft magnetic 
photons. A quantum kinetic equation for hard fermions is obtained directly in 
real time from a field theoretical approach improved by the dynamical renor- 
malization group. The collision kernel is time- dependent and infrared finite. 
In linear response the fermion distribution function relaxes with an anomalous 
exponential law with an exponent twice as large as that for the mean field. 
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I. INTRODUCTION 



The study of nonequilibrium phenomena under extreme conditions play a fundamental 
role in the understanding of ultrarelativistic heavy ion collisions and early universe cosmol- 
ogy. Forthcoming relativistic heavy ion experiments at BNL Relativistic Heavy Ion Collider 
(RHIC) and CERN Large Hadron Collider (LHC) aim to search for a deconfmed phase 
of quarks and gluons, the quark gluon plasma (QGP), which is predicted by lattice QCD 
simulations [ffLQ to emerge at a temperature scale T > 200 MeV. 

Recent results from CERN Super Proton Synchrotron (SPS) seem to confirm the 
main theoretical ideas that in the central region in ultrarelativistic heavy ion collision a de- 
confined plasma of quarks and gluons forms which expands and cools rapidly and eventually 
hadronizes. 

Current estimates based on energy deposited in the central region for -y/i / nucleon-pair ~ 
200 GeV at RHIC suggest that the lifetime of a deconfmed QGP is of order 10-50 fm/c 
At such unprecedented short time scales, an important aspect is an assessment of thermal- 
ization time scales and the potential for non-equilibrium effects associated with the rapid 
expansion and finite lifetime of the plasma and their impact on experimental observables. 
Lattice QCD is simply unable to deal with these questions because simulations are restricted 
to thermodynamic equilibrium quantities and a field-theoretical nonequilibrium approach is 
needed for an accurate description of the formation and evolution of the QGP. An important 
and pioneering step in this direction was undertaken by Geiger who applied transport 
methods combined with perturbative QCD (pQCD) to obtain a quantitative picture of the 
evolution of partons in the early stages of formation and evolution of the plasma. 

The consistent study of the evolution of partons in terms of pQCD cross sections that 
include screening corrections to avoid the infrared divergences associated with small angle 
scattering lead to the conclusion that quarks and gluons thermalize on time scales of a 
few fm/c ||. The necessity of a deeper understanding of equilibrium and non-equilibrium 
aspects of the quark-gluon plasma motivated an intense study of the abelian and non-abelian 
plasmas in extreme environments. A major step towards a consistent description of non- 
perturbative aspects was taken by Braaten and Pisarski who introduced a novel 



resummation method that re-organizes the perturbative expansion in terms of the degrees 
of freedom associated with collective modes Q, rather than bare particles. This program, 
called the hard thermal loop or HTL program is now at the heart of most treatments of 
equilibrium aspects of abelian and non-abelian plasmas ]nfl . 

Thermal field theory provides the tools to study the properties of plasmas in equilib- 
rium [ Tl| — p~4 1 , but the consistent study of nonequilibrium phenomena in real time requires 



the methods of nonequilibrium field theory |TJJ (see also Refs. |n],|T7| f° r further references). 
The study of the equilibrium and nonequilibrium properties of abelian and non-abelian 
plasmas as applied to the QGP has as ultimate goal a deeper understanding of the poten- 
tial experimental signatures of the formation and evolution of the QGP in ultrarelativistic 
heavy ion collisions. Amongst these, photons and dileptons (electron and/or muon pairs) 
produced during the early stages of the QGP are considered as some of the most promis- 
ing sig nals PHU]. Since photons and lepton pairs interact electromagnetically their mean 
free paths are longer than the estimated size of the QGP fireball ~ 10 — 50 fm and un- 
like hadronic signals they do not undergo final state interactions. Therefore photons and 
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dileptons produced during the early stages of QGP carry clean information from this phase. 

The goals of this work. In this work we aim to provide a comprehensive study of 
several relevant aspects of the nonequilibrium dynamics of an abelian QED plasma in real 
time. Many features of full QCD are similar to those of the abelian (QED) theory P, IT], 9] , 



in particular leading contributions in the HTL limit can be straighfordarly generalized from 
QED to QCD, thus the leading results for relaxation and photon production from a QGP 
|19| , |20[| can be understood from the study of a QED plasma. 

We utilize a gauge invariant formulation, available in abelian gauge theories that cir- 
cumvents the possible ambiguities associated with gauge invariance |H],[np . Thus gauge and 



fermionic mean fields and distribution functions are automatically gauge invariant. We im- 
plement and apply the method of the dynamical renormalization group, introduced recently 
to study non-equilibrium phenomena directly in real time l2l|j22|j to extract a consistent, 
non-perturbative description of real-time dynamics out of equilibrium. 
In particular we focus on the following: 

(i) The real time evolution of gauge mean fields in linear response in the HTL approxima- 
tion. The goal here is to study directly in real time the relaxation of (coherent) gauge 
field configurations in the linearized approximation to leading order in the HTL pro- 
gram. Whereas a similar study has been carried out in scalar QED |[23] and confirmed 



numerically in |24|] the most relevant case of spinor QED has not yet been studied in 
detail. 

The quantum kinetic equation that describes the evolution of the distribution function 
of photons in the medium, again to leading order in the HTL approximation. This 
aspect is relevant to study photon production via off-shell effects directly in real time. 
As explained in detail, this quantum kinetic equation, obtained from a microscopic field 
theoretical approach based on the dynamical renormalization group pT| , p2| displays 
novel off-shell effects that cannot be captured via the usual kinetic description that 
assumes completed collisions |18^p0| . 



(iii) The evolution in real time of fermionic mean fields features anomalous relaxation aris- 
ing from the emission and absorption of magnetic photons (gluons) which are only 



dynamically screened by Landau damping |10| , P5| , ^6| . The fermion propagator was 
studied previously in real time in the Bloch-Nordsieck approximation which provides 
a resummation of the infrared divergences associated with soft photon (or gluon) 



bremsstrahlung in the medium pB] , pB| . In this article we implement the dynamical 
renormalization group to study the evolution of fermionic mean fields providing an 
alternative to the Bloch-Nordsieck treatment. 

(iv) We obtain the quantum kinetic equation for the fermionic distribution function for 
hard fermions via the implementation of the dynamical renormalization group. There 
has recently been an important effort in trying to obtain the effective kinetic (Boltz- 



mann) equations for hard charged (quasi)particles p7|-p9| but the collision kernel in 
this equation features the logarithmic divergences associated with the emission and 
absorption of soft magnetic photons (or gluons) p{|. The dynamical renormaliza- 



tion group leads to a quantum kinetic equation directly in real time bypassing the 
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assumption of completed collisions and leads to a time- dependent collision kernel free 
of infrared divergences. 



Summary of the main results. The main results of this study are summarized as 
follows. 

Relaxation of gauge mean fields. We studied the relaxation of a gauge mean field 
in linear response to leading order in the HTL approximation both for soft momentum 
k < eT and for semihard momentum eT <C k <C T under the assumption of weak 
electromagnetic coupling. 

Soft momentum (k ~ eT): in this case the relaxation of the gauge mean field is dom- 
inated by the end-point contribution of the Landau damping cut. As a consequence, 
the soft gauge mean field relaxes with a power law long time tail of the form 



a T (k,t) fe ^ >1 a r (k,0) 



k 2 Z T (k) 
uj^{k) 



cos[u>T(k)t] 



12 cos kt 
e 2 T 2 t 2 



where ur(k) is the transverse photon pole and Z^{k) is the corresponding residue. 
We note that in spite of the power law tail the gauge mean field relaxes towards the 
oscillatory mode determined by the transverse photon pole. This reveals that the soft 
collective excitation in a plasma is stable in the HTL approximation. 

Ultrasoft momentum (k <C eT): In the region of ultrasoft momentum the spectral den- 
sity divided by the frequency features a sharp Breit-Wigner peak near zero frequency 
in the region of Landau damping, with width = 12/c 3 /7re 2 T 2 . We find that the am- 
plitude of a mean field of transverse photons prepared via a source that is adiabatically 
switched-on is given by 



, fct>i,r\t<i 
a r (M) = ~ a T (k,0) 



k 2 Z T {k) 



cos[cj T (A;)t] + e~ Tkt 



We emphasize that the exponential decay is a consequence of the sharp resonance near 
zero frequency in the Landau damping region of the spectral density and only arises 
if the mean field is prepared by an external source whose time Fourier transform has 
a simple pole at zero frequency, such is the case for an adiabatically prepared initial 
state. This result confirms those found numerically in Ref. [^4 . 

Semi-hard momentum (eT <C k -C T): In this region both the HTL approximation 
and the perturbative expansion are formally valid. However the spectral density in 
the Landau damping region is sharply peaked near u = k and the transverse photon 
pole approaches the edge of the Landau damping region from above. We find that 
although the perturbative expansion is in principle valid, the sharp spectral density 
near the edge of the continuum results in a breakdown of the perturbative expansion. 
The dynamical renormalization group provides a consistent resummation of the lowest 
order HTL perturbative contributions in real time, leading to the following relaxation 
of the mean field at intermediate asymptotic times: 



a T (k,t) 



fct>i 



12k' 2 
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where tq ~ 1/k and «4.T(k, t) is an oscillating function. The anomalous exponent is a 
consequence of an infrared enhancement arising from the sharp spectral density near 
the threshold of the Landau damping region for semihard momentum. The crossover 
to exponential relaxation due to collisional processes at higher orders is discussed. 

Quantum kinetic equation for the photon distribution function. Using the 
techniques of nonequilibrium field theory and the dynamical renormalization group, we 
obtain the quantum kinetic equation for the distribution function of semihard photons 
eT k ^ T to lowest order in the HTL approximation assuming that the fermions 
are thermalized. An important result is that the collision kernel is time- dependent and 
the dynamical renormalization group reveals that detailed balance emerges during 
microscopic time scales, i.e, much shorter than the relaxation scales. In the linearized 
approximation we find that the departure from equilibrium of the photon distribution 
function relaxes as: 

e 2 T 2 

5nl(t) = Snl(t ) (-^) ^ for k(t - t ) > 1 , 

where To ~ and to is the initial time. Furthermore, this quantum kinetic equation 
allows us to study photon production by off-shell effects, which to leading order in 
the HTL approximation are determined by photon bremsstrahlung and is of order 
a. Extrapolating the result from QED to thermalized QGP with two flavors of light 
quarks, we find that the total number of hard photons at time t per invariant phase 
space volume to lowest order is 

5aT 3 

^) = 7^^{ ln W-*o)]+7£-l} for fc(*-*o >1- 

with to ~ lfm/c is the time scale at which the QGP plasma is thermalized. 

We find that for a quark-gluon plasma at temperature T ~ 200 MeV and of life- 
time 10 < (t — t Q ) < 50 fm/c, the hard (k ~ T) photon production by off-shell 
bremsstrahlung is comparable to that from Compton scattering and pair annihilation 
of order a a s [0,^]] • 

Relaxation of fermion mean fields. We implement the dynamical renormaliza- 
tion group resummation to study the real-time relaxation of a fermion mean field for 
hard momentum. The emission and absorption of magnetic photons which are only 
dynamically screened by Landau damping introduce a logarithmic divergence in the 
spectral density near the fermion mass shell. The dynamical renormalization group 
resums these divergences in real time and leads to a relaxation of the fermion mean 
field for hard momentum given by: 

</>(k,t) e -^[ln M+ 0.12652...] x osciUating phaseS; 

with up being the plasma frequency. 

Quantum kinetics for the fermion distribution function. We obtain a quantum 
kinetic equation for the distribution function of hard fermions using non-equilibrium 
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field theory and the dynamical renormalization group resummation. Assuming ther- 
malized photons we find that the collision kernel is infrared finite but time- dependent. 
In the linearized approximation the distribution function relaxes as: 

5ni(t) kt ^ 5ni(t ) e -2^-*o)[M^)+o.i2652...] ? 
where the anomalous relaxation exponent is twice that of the mean field. 

The article is organized as follows. In Sec. II we review the main ingredients of nonequi- 
librium field theory, the initial value problem formulation for relaxation of nonequilibrium 
mean fields, and the dynamical renormalization group approach to quantum kinetics. In 
Sec. Ill we first study relaxation of the photon mean field in the hard thermal loop limit for 
soft k < eT and semihard eT <C k <T photon momentum and obtain the quantum kinetic 
equation for the (hard and semihard) photon distribution function. In Sec. IV we study re- 
laxation of fermionic mean fields for hard momentum and the quantum kinetic equation for 
the fermion distribution function. Our conclusions and some further questions are presented 
in Sec. V. 



II. GENERAL ASPECTS 

As stated in the introduction our goal is to provide a systematic study of non-equilibrium 
phenomena in a hot abelian plasma. In particular we focus on a detailed study of the real 
time relaxation of expectation values of fermions and gauge bosons, i.e, fermionic and gauge 
mean fields as well as the quantum kinetics for the evolution of the expectation value of 
the number operator associated with fermions and photons. A fundamental issue that must 
be addressed prior to setting up our study is that of gauge invariance. In the abelian 
case it is straightforward to reduce the Hilbert space to the gauge invariant subspace and 
to define gauge invariant charged fermionic operators. The description in terms of gauge 
invariant states and operators is best achieved within the canonical formulation which begins 
with the identification of the canonical field and conjugate momenta and the primary and 
secondary first class constraints associated with gauge invariance. The physical states are 
those annihilated by the constraints and physical operators commute with the (first class) 
constraints. This program has been implemented explicitly in the case of scalar quantum 



electrodynamics p2| , |30| and the fermionic case can be treated in the same manner with few 
minor technical modifications. 

The final result of this formulation is that the Hamiltonian acting on the gauge invariant 
states and written in terms of gauge invariant fields is exactly equivalent to that obtained in 
Coulomb gauge, which is the statement that Coulomb gauge describes the theory in terms 
of the physical degrees of freedom. Furthermore the instantaneous Coulomb interaction can 
be traded for a Lagrange multiplier [52] leading to the following Lagrangian density 



£ = ^ {ifi - e 7o A + e 7 ■ A T - m) + - [(<^A T ) 2 + (VA ) 2 

where A^ is the transverse component of the vector potential and Aq is not to be interpreted 
as the time component of the gauge vector potential but is the Lagrange multiplier associated 
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with the instantaneous Coulomb interaction. We emphasize that the fermionic charged fields 
\I/ as well as the Lagrange multiplier A are gauge invariant fields (see Refs. [ 22 J30| ) . 



The fully renormalized equations of motion for the mean fields can be obtained by fol- 



lowing the formulation presented in Refs. |31] , |32| . However, the counterterms that eliminate 



the zero temperature ultraviolet divergences are temperature independent and play 
no role in the present context, thus we will neglect the zero temperature renormalization 
counterterms in our study. 

Furthermore we consider a neutral plasma (i.e, with zero chemical potential for the 
charged fields) at a temperature T ^> m with m the (renormalized) mass of the fermions, 
hence in what follows we neglect the fermion mass unless otherwise stated. 

A. Nonequilibrium Field Theory 

The formulation of nonequilibrium quantum field theory in terms of the Schwinger- 
Keldysh or the closed-time-path (CTP) is standard UlSHT^I. A path integral representation 
requires a contour in the complex time plane, running forward and then backwards in time 
corresponding to the unitary time evolution of an initially prepared density matrix, with an 
effective Lagrangian in terms of fields on the different branches of the path 



£ noncq = tf+ A+, A+] Aj, A 



o . 



where the "+" ( "— " ) superscripts for the fields refer to fields defined in the forward (back- 
ward) time branches. This description leads to a straightforward diagrammatic expansion 
of the nonequilibrium Green's functions in terms of real-time propagators but with modified 
Feynman rules (as compared to standard field theory) [[15-17]. The free fermion and photon 
propagators are given in detail in the Appendix. 

Since we study two different type of situations: (i) the relaxation of expectation values 
in the linearized approximation (linear response) and (ii) the kinetic equation that describes 
the evolution of a non-equilibrium distribution function, we make a distinction between the 
real-time propagators used in each case. 

(i) In the case of relaxation of expectation values in linear response, the distribution 
functions that enter in the fermion and photon propagators (|A1|)- (|A"3D are those with 
equilibrium distributions. 

(ii) When we study the kinetic equations for the expectation value of the fermion number 
and photon number operators, we assume that the initial density matrix is diagonal in 
the basis of the quasiparticles (see below) but with nonequilibrium initial distributions. 
Thus the propagators obtained with this initial density matrix have the free-field form, 
as given in (|A1|) - (|A^) but the initial occupation numbers are not the equilibrium Fermi- 
Dirac or Bose-Einstein distribution. 



B. Real-time relaxation of mean fields in linear response 

The mean fields under consideration are the expectation value of either fermion or gauge 
field operators in the nonequilibrium state induced by the external source. Obviously in 
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equilibrium and in the absence of external sources these must vanish, and we introduce 
time-dependent sources to induce an expectation value of these operators. Our strategy 
to study the relaxation of these mean fields as an initial value problem is to prepare these 
mean fields via the adiabatic switching-on of an external current f3lH33f. Once the current 
is switched off the expectation values must relax towards equilibrium and we study this 
real-time evolution. This formulation to study the real-time evolution of mean fields has the 
advantage that it leads to a straightforward and systematic implementation of the dynamical 
renormalization group method as explained in detail in Refs. 

The initial value problem formulation for studying the linear relaxation of the mean fields 
begins by introducing c-number external sources coupled to the quantum fields. Let r)(x) 
be the Grassmann- valued fermionic source and Jt{x) be the electromagnetic sources,[] then 
the Lagrangian density becomes 

£ -> C + vi/77 + r/^ + J T • A T . 

The presence of external sources will induce responses of the system. Let Q(x) be a generic 
quantum field (i.e., fermionic or bosonic) with $ ± (x) the fields defined on the forward and 
backward branches, respectively, and J{x) the corresponding c-number external source (the 
same for both branches). The expectation value of $(x) induced by J(x) in a linear response 
analysis is given by 



(x,t) = ($ ± (x,t)) J 

r r+oo 

= / dV / dt' G rct (x, t; x', if) J(x', if) 



with the retarded Green's function 

G re t(x, t- x', t')=i\ ($+(x, i)$ + (x', t')) - ($+(x, t)<i-(x', If] 



$(x,t),$(x',i')J T /^-0, 

where expectation values are computed in the CTP formulation in the full interacting theory 
but with vanishing external source, 6(t — t') is the step function and the subscripts =F refer 
to the commutator (— ) for bosonic or anticommutator (+) for fermionic fields. 

A practically useful initial value problem formulation for the real-time relaxation of mean 
fields is obtained by considering that the external source is adiabatically switched on in time 
at t = — oo and suddenly switched off at t = t , i.e., 

J(x, t) = J(x) e e( *-* o) 0(t o - t) , e -> 0+ . (2.1) 

The adiabatic switching on of the external source induces a mean field that is dressed 
adiabatically by the interaction. Then for t > to, after the external current has been 



x In this article we will not discuss relaxation of the longitudinal gauge field Ao(x) associated with 
the instantaneous Coulomb interaction, hence the corresponding external source is neglected. 



S 



switched off, the mean field will relax towards equilibrium, and our aim is to study this 



relaxation directly in real time |33 



The equations of motion of the mean fields are obtained via the tadpole method |T7] and 



are automatically causal and retarded. The central idea of the tadpole method is to write 
the field into a c-numbered expectation value plus fluctuations around it, i.e., writing 

$ ± (x,t) = 0(x,t) + X ± (x,t) with ($ ± (x,t)) J = 0(x,t) . 

The equation of motion for 0(x, t) is obtained by requiring (x ± (x, t)) = to all orders in 
perturbation theory. 

For the study of relaxation of the mean fields, we assume that at t = — oo the system 
is in thermal equilibrium at a temperature T and henceforth choose t Q = for convenience. 
The retarded and the equilibrium (i.e., time translational invariant) nature of G re t(x, t; x', t') 
and the adiabatic switching on of J(x, t) entail that 



(x,t = 0) =0 o (x) 
(x,t < 0) = , 



(2.2) 



where 0o( x ) is determined by J(x) (or vice versa). It is worth pointing out that <j>(x,t) 
is not specified at t = although 0(x, t < 0) = 0. This is because the external source is 
switched off at t = 0.0 



C. Main ingredients for quantum kinetics 



In Ref. |22| a systematic treatment to obtain the kinetic equation that determines the 
real-time evolution of the distribution function (defined as the expectation value of the num- 
ber operator) was established, leading to its interpretation as a dynamical renormalization 
group equation. Here we summarize the basic steps to obtain the relevant quantum kinetic 



equations from first principles and we refer the reader to Ref. |22| for more details. 



(i) Identify the proper degrees of freedom (quasiparticles) to be described by the kinetic 
equation, the corresponding microscopic time scale r micro associated with their oscil- 
lation, and the number operator N(k, t) that counts these degrees of freedom with 
momentum k at time t. 

ii) Use the Heisenberg equations of motion to derive a general rate equation for dn-k(t)/dt 
with rik(t) = (N(k, £)), where the expectation value is taken in the initial density ma- 
trix. Assuming that the initial density matrix is diagonal in the basis of free quasipar- 
ticles but with nonequilibrium distribution functions, we can expand the rate equation 
in perturbation theory using nonequilibrium Feynman rules and free real-time propa- 
gators in terms of nonequilibrium distribution functions. Since the resulting expression 



2 There could be initial time singularities associated with our choice of the external source (see 
Ref. []33| ) . However the long time asymptotics is not sensitive to the initial time singularities |33| , 
and we will not address this issue here as it is not relevant for the asymptotic dynamics. 
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is a functional of the distribution function at initial time, the solution of the kinetic 
equation can be obtained by direct integration. 

This solution is typically characterized by the emergence of secular terms, i.e, terms 
that grow in time in the intermediate asymptotic regime (r micro <C t <C r re i) where the 
scale r re i is that at which perturbation theory breaks down |2T]f22|] . 

(iii) These secular terms growing in time become divergent if the perturbative solution is 
extrapolated to long times. However, perturbation theory is valid and dominated by 
the secular terms during the intermediate asymptotic time scales r micro <C t <C r re i. 
The dynamical renormalization group is used to resum the secular terms through a 
renormalization of the initial distribution functions. This introduces an arbitrary time 
scale t which serves as a renormalization point. The r-independence of the renormal- 
ized solution leads to the dynamical renormalization equation |21| , |22| . This dynamical 
renormalization equation describing the evolution of the quasiparticle distribution is 
recognized as the quantum kinetic equation. 



III. PHOTONS OUT OF EQUILIBRIUM 

In this section we study non-equilibrium aspects of photon relaxation and production. 
We begin by analyzing the relaxation in real time of a photon condensate in linear response 
both in the case of soft k < eT and semihard (or semisoft) eT <C fc C T assuming the 
electromagnetic coupling to be small. We then continue with a study of the production of 
semihard photons via off-shell effects. 



A. Relaxation of the gauge mean field 

We begin with the relaxation of soft photons of momenta k ~ eT. As mentioned in 
the previous section, the equation of motion for the transverse photon mean field can be 
derived from the tadpole method [17j by decomposing the full quantum fields into c-number 



expectation values and quantum fluctuations around the expectation values: 

A^ (x, t) = a T (x, t) + A£ (x, t) , 
a T (x,t) = (A±(x,f)> . 

With the external source Jt( x , t) of the adiabatic form ( |2.1| ) chosen to vanish at t > 0, we 
find the linear equation of motion for the spatial Fourier transform of ay(x, t) for t > to 
be given by 



where 



(J^ + fc2 ) a ^( k > + f dt ' n ^( k > * - M k > = 

a T (k,t) = J d 3 x e" ik x a T (x,t) . 



(3.1) 



10 



Here ILr(k, t — t') is the transverse part of the retarded photon self-energy and we have 
made explicit use of its retarded nature in writing Eq. ( [3.1|) , which evidently is causal and 
retarded. 

Using the nonequilibrium Feynman rules and the real-time free fermion propagators given 
in the Appendix, we find the transverse polarization to be given by 

/+oo 
du p T (uj, k) sin[o;(t - t')} , (3.2) 
-oo 

in terms of the spectral density to one-loop order 

pr(u, k) = -e 2 J J^L{[<f(u; ~ V ~ q) ~ 5(u + p + q)} [1 + (k • p)(k ■ q)] [1 - n F (p) - n F {q)\ 
+ [6(u -p + q)-5(u+p-q)} [n F (q) - n F (p)] [1 - (k ■ p)(£ ■ q)]} , (3.3) 
where p = k + q. 

It is now convenient to introduce an auxiliary quantity 7r^(k, t — £') as follows 

U T (k,t-t') = ^n T (k,t-t'). 

and hence, 

/+OO fjjjj 
— pr(u, k) cos[uj(t - t')] . 
-oo CO 

Upon integration by parts the equation of motion takes a form that displays clearly the 
nature of the initial value problem: 

+ k 2 ^j a T (k, t) - jT* dt' vr T (k, t - t') a T (k, t') + vr r (k, 0) a T (k, t) = , (3.4) 

where a dot denotes derivative with respect to time and we used that a^(k, t) — for t < 0. 
Eq. (3.4) together with the initial conditions specified at t — yields a well-defined initial 
value problem. As discussed above, the initial value a T (k, t = 0) = a T (k, 0) is determined 
by Xr(k, i), and we choose 

a T (k,t = 0)=0 (3.5) 



consistently with the adiabatic switching-on condition (|2.2j ). 

Eq. (3.4) can be solved by the Laplace transform method. We define, 

POO 

a T (s, k) = / dt e~ st a T (k,t) , 
Jo 

roo 

7Pr(s,k)=/ dt e" s * 7r r (k,t) , (3.6) 
Jo 

with Re s > 0. In terms of which the Laplace transform of the equation of motion leads to 
s 2 + k 2 + ft r (s, k)l a r (s, k) = [s - 7f T (s, k)] a r (k, 0) , (3.7) 
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where ILr(s, k) is the Laplace transform of Hr(k, t — t')\ 

n T (s,k) =vr T (k,0)- S 7f T (s,k) . 
The solution of Eq. ( |3.7|) is given by 

a r (a, k) = i {l - A T (s, k) [k 2 + vr T (k, 0)]} a T (k, 0) , (3.8) 

where the retarded transverse photon propagator Ay(s, k) is given by 

A r (s,k) = [s 2 + A; 2 + ri T ( S ,k)]- 1 . (3.9) 

The real-time evolution of a^(k, t) is obtained by performing the inverse Laplace transform 
along the Bromwich contour which is to the right of all singularities of a^(s, k) in the complex 
s-plane. 

We note that this result is rather different from that obtained in Ref. [2iJ in the structure 
of the solution, in particular the prefactor 1/s in Eq. ( |3.8|) . This can be traced back to the 
different manner in which we have set up the initial value problem, as compared to the case 
studied in [p3] . The adiabatic source ( |2.1| ) (for £ = 0) has a Fourier transform that has 
a simple pole in the frequency plane, this translates into the 1/s factor in Eq. ( |3.8| ). This 
particular case was not contemplated in those studied in Ref. [^3] since the source ( j2.1Q is 
not a regular function of the frequency. This difference will be seen to be at the heart of 
important aspects of relaxation of the mean field condensate in the soft momentum limit as 
discussed in detail below. 

It is straightforward to see that the residue vanishes at s = in the solution (|3.8|) . The 
singularities of a^(s, k) are those arising from the retarded transverse photon propagator 
A T (s,k). 



1. Soft photons k ~ eT: real-time Landau damping 



For soft photons of momenta k ~ eT, the leading (~ e 2 T 2 ) contribution to the pho- 
ton self-energy arises from loop momenta of order T and is referred to as a hard thermal 
loop (HTL) [§,0. In this region of momenta, the contribution of the fermion loop is non- 
perturbative. 

In the HTL approximation (s ~ k <C q), after some algebra we find that the HTL 
contribution to Il r (s, k) arises exclusively from the terms associated with 5(u> =F V =t q) m 
Eq. (|3.3|) and reads 



e 2 T 2 
~12 



n r (s,k) 



k 2 j is — k k 2 



(3.10) 



The delta functions S(u =F p ± q) have support below the light cone (i.e., uo 2 < k 2 ) and 
correspond to Landau damping processes |J in which the soft photon scatters a hard fermion 
in the plasma. 

The analytic continuation of Ht{s, k) in the complex s-plane is defined by 
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n T (cu, k) = n T (s = -iu ± + , k) = Ren T (cu, k) ± i Imn T (c;, k). 



(3.11) 



From Eq. fl3.2|) it is straightforward to see that the real and imaginary parts are related by 
the usual dispersion relation. The analytical continuation of A T (s, k) and 7fr(s,k) can be 
defined analogously, and they are related to n T (u;, k) by 

A t (uj, k) = -[u 2 - k 2 - U t {uj, k)]" 1 , 
ReHr(u;, k) = ^(k, 0) — uj Ini7rr(u;, k) , 
ImHr(u;, k) = uj Re7Pr(co>, k) . 

The analytical continuation of ITt(s, k) in the HTL limit thus reads 



n T (cu, k) 



e 2 T 2 



12 



w UJ 

~k? + k 



UJ 



In 



uj + k 



— i- 



ne 2 T 2 uj 
12 k 



to 



uj — k 
uj 2 ). 



In the HTL limit it is straightforward to prove that vrr(k, 0) = 0. We note that ImHr(co>, k) 
only has support below the light cone and vanishes linearly as uj — > k from below. Since for 
soft momenta k -C eT the HTL photon self-energy is comparable in magnitude to or larger 
than the free inverse propagator, it has to be treated nonperturbatively. 

Isolated poles of At(w, k) in the complex tu-plane correspond to collective excitations. 
The transverse photon pole ujT{k) is real and determined by |[TT| 

uj 2 T {k) - k 2 - ReIl T [uj T (k), k] = . 



For ultrasoft photons k -C eT, the dispersion relation of the collective excitations is |TT 

fc 4 \ 



uj T (k) 



6- 
-i 

5 



oj 2 p +-k 2 + 



e 2 T 2 



where ujp = eT/3 is the plasma frequency. Consequently, the retarded transverse photon 
propagator Ay(s, k) has two isolated poles at s = ±iujT{k) and a branch cut from s = —ik 
to s = ik. 

Having analyzed the analytic structure of At{s, k), we can now perform the integral along 
the Bromwich contour to obtain the real-time evolution of a-r(k, t). Closing the contour in 
the half-plane Res < 0, we find ar(k, t) for t > to be given by 



a T (k,t) = a^ e (M + a^(k,t) 



with 



pole 



(M) 



k 2 Z T (k) 
uj^(k) 



cos[ujT(k)t] a r (k, 0) 



— (3 T {uj,k) e- wt a T (k,0), 

-k UJ 



(3.12) 
(3.13) 



where 
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Z T (k) 
p T (u,k) 



1 - 



dReU T (cu, k) 
dcu 2 



: 2 T 2 
12 



k V fc 2 j 



A; 2 



g 2 T 2 
12 



12 " fc V 1 k 2 



(3.14) 



The solution evaluated at t = must match the initial condition, this requirement leads 



immediately to the sum rule [11 



+°°dq 1 
-oo g <r 



where prilo, q) is the HTL-resummed spectral density for the transverse photon propagator :[] 



Pt(Qo,q) = -ImA T (g , g) 



7T 



sgn(g ) Z T (g) 5[g 2 - 4(g)] + /3 T (g , g) 9(q 2 - g 2 ) 



(3.15) 



A noteworthy feature of the cut contribution ( [3.13 ) is the factor u in the denominator. The 
presence of this factor can be traced back to the preparation of the initial state via the 
adiabatically switched-on external source, which is the switched off at t = 0. The Fourier 
transform of this source is proportional to 1/w and results in the prefactor of /3t(uj, k) in 
Eq. (EEII). 

For k ~ eT, Eq. ( |3.13| ) cannot be evaluated in closed form but its long time asymptotics 
can be extracted by writing the integral along the cut as a contour integral in the complex 
k>plane. The integration contour C is chosen to run clockwise along the segment —k<uo<k 
on the real axis, the line u = k — iz with < z < oo, then around an arc at infinity and 
back to the real axis along the line u = —k — iz. After some algebra we obtain the following 
expression 



SL^ ut {k,t) = k 2 le~ ikt J o dze 



ifa^ — iz, k) 



— 2ni Res 



poles 
inside C 




Because of the exponential factor e~ zt in the integrals , the dominant contributions to the 
integral at long times (t ^> 1/k) arise from the end points of the branch cut with z <C k, 
and lead to a long time behavior characterized by a power law: 



fc*>i 



12 cos kt 
e 2 T 2 t 2 



a T (k, 0) 



(3.16) 



3 In our notation the spectral density for the self-energy is denoted by p, and the spectral density 
for the corresponding propagator is denoted by p. 
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which confirms the results of Ref. ||23|| . The contribution from the poles inside the contour 



C is dominated at long times by the pole closest to the real axis. 

For very soft momenta k <C eT the function /3t(o;, k)/uo is strongly peaked at u = as 
depicted in Fig. |l|. For k -C eT and k ^> u, Eq. (|3.14j) takes a Breit-Wigner form 

-f3 T {uj } k) = with T fc = ^— . 3.17 

Using this narrow-width approximation to evaluate Eq. ( |3.13| ) yields, 

a r (k,0). (3.18) 



i^(k,t) kt ^ 



The end-point contribution which results in a power law as in Eq. fl3.16|) is very small com 



pared with Eq. fl3.18|) for time t < l/T^. This power law becomes dominant at time scales 
much longer than l/T^, because its amplitude at t ~ l/r^ is of order (k/eT) 6 -C 1 in the 
soft momentum limit. 

Thus, for adiabatically prepared mean fields of soft momentum fc < eT we find the 
real-time behavior 

i fct>i,r fc t<i , (k 2 Z T (k) , ,,. , r .1 

a r (k,t) = ~ a T (k, 0) | ; cosMfc)t] + e- rfct | . 

It is worth emphasizing that the exponential decay is not the same as the usual decay 
of an unstable particle (or even collisional broadening) for which the amplitude relaxes to 
zero at times longer than the relaxation time. In this case the asymptotic behavior of the 
amplitude is completely determined by the transverse photon pole u>T{k), and to this order 
in the HTL approximation the collective excitation is stable. The exponential relaxation 
does not arise from a resonance at the position of the transverse photon pole but at zero 
frequency. Clearly we expect a true exponential damping of the collective excitation at 
higher order result of collisional broadening;. 

Our results thus confirm those found in hot scalar quantum electrodynamics |23] wherein 



a numerical analysis of the relaxation of the mean field was carried out. In that reference the 
authors studied the real time evolution in terms of local Hamiltonian equations obtained by 
introducing a non-local auxiliary field. The initial conditions chosen there for this non-local 
auxiliary field correspond precisely to the choice of an adiabatically switched on current 
leading to an adiabatically prepared initial value problem. 

We emphasize that the exponential relaxation ( |3.18| ) can only be probed by adiabatically 
switching on an external source [see Eq. ( |2.1|) 1 whose temporal Fourier transform has a 
simple pole at ui = which is the origin of the prefactor 1/cu in Eq. (|3.13| ). That is the 
adiabatic preparation of the initial state excites this pole in the Landau damping region. 
For external sources whose temporal Fourier transform is regular at u = no exponential 
relaxation arises, in agreement with the results of Ref. . Thus we find that the exponential 
relaxation is not a generic feature of the evolution of the mean field, but emerges only for 
particular (albeit physically motivated) initial conditions. 
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2. Semi-hard photons eT <C k <C T: anomalous relaxation 



Most of the studies of the photon self-energy in the hard thermal loop approximation 
focused on the soft external momentum region k <C eT (with e C 1) |6|- [Tl|1 wherein the 
contribution of the hard thermal fermion loop must be treated non-nonperturbatively. 

However there are important reasons that warrant a consideration of the region of semi- 
hard photons eT < C T. From a phenomenological standpoint, hard and semihard 
photons produced in the QGP phase are important electromagnetic probes of the decon- 
fined phase |l|J^, p^ , ^0[1 hence the study of all aspects of the space-time evolution, production 
and relaxation of photons to explore potential experimental signatures is warranted. 

A more theoretical justification of the relevance of this region is that whereas the photon 
self-energy for this region of momenta is still dominated by the hard thermal fermion loop 
contribution ( |3.10| ), its contribution to the full photon propagator is now perturbative as 



compared with the free inverse propagator. The validity of perturbation theory in this 



regime licenses us to study the real-time evolution with the tools developed in refs. pl|j22 
that provide a consistent implementation of the renormalization group to study real-time 
phenomena. 

The dynamical renormalization group approach introduced in Refs. pl] , ^2| is particularly 
suited to study the real-time evolution in the case in which there are (infrared) threshold 
singularities in the spectral density. To understand the potential emergence of anomalous 
thresholds in the semihard and hard limit for the photon self-energy we note that at leading 
order in HTL the Laplace transform of the inverse photon propagator Eqs. ( |3 . 9| ) - (|37T0| ) reads 

e 2 T 2 s 2 



A^fok) = (s 2 + k 2 



e 2 T 2 is is + k 

1 H -r^- In - 

12 k 3 is-k 



6 k 2 



The transverse photon poles are for s = ±i[k + e 2 T 2 /l2k + 0(T /k 3 )] which in the limit 
k ^ eT approach s — > ±ik. That is, 



UT (k) = k + — - + O — . (3.19) 




In this limit the poles merge with the thresholds of the logarithmic branch cut from Landau 
damping and are no longer isolated from the continuum. The approach of the pole to the 
Landau damping threshold and the enhancement of the spectral density near threshold for 
k ^> eT is clearly displayed in Fig. |2|. While the perturbative expansion in the effective 
coupling e 2 T 2 /k 2 is warranted in the semihard case under consideration, the wave function 
renormalization constant evaluated as the residue at the pole is a non-analytic function of 
this coupling and given by 



e 2 T 2 



k 24? +3 + ° h^T. 



This is obviously a consequence of the logarithmic threshold singularities in the semihard 
regime. This threshold singularity is reminiscent of those studied in detail in Ref. [^TJ. In 



this reference it was shown that the Bloch-Nordsieck resummation which is equivalent to a 
renormalization group improvement of the space-time Fourier transform of the self-energy, 
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leads to a real-time evolution which is obtained by the implementation of the dynamical 

Since the infrared threshold singularities in the 
we now implement the dynamical 



renormalization group in real time pi 
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semihard limit are akin to those studied in Ref. 
renormalization group program advocated in that reference to study the relaxation of the 
photon mean field in the semihard limit directly in real time. 

Perturbation theory in terms of the effective coupling e 2 T 2 /k 2 is in principle reliable for 
semihard photons, hence we can try to solve Eq. ( ft.4|) by perturbative expansion in powers 
of e 2 (the true dimensionless coupling is e 2 T 2 /k 2 ). Writing 



a T (M) 
tt t (M-0 



a5? ) (k,t) + e 2 a^ 1) (k,t) + C(e 4 ) 
e 2 Tr£\k,t-t')+0(e 4 ), 



and expanding consistently in powers of e 2 we obtain a hierarchy of equations: 
( d 2 



dt 2 



+ k 2 1 - (0) 



a^ ; (k,t) = 0, 



4 1] (k, t) = f dt' 4 1} (k, t - t') 4 0) (k, , 
Jo 



where we have explicitly used the fact that 7r^(k, 0) = in the HTL limit. Starting from 
the solution to the zeroth-order equation 



4 0) (k,t) = J2 K(k) e~ ikt + A*(k) e ikt ] £ A (k) 



A=l 



(3.20) 



with £ A (k) being the polarization vector, and the retarded Green's function of the unper- 
turbed problem: 

GS(M-*') = Sin[ t~ t ' )1 Ht -if), 



the solution to the first-order equation reads 

2 r+oo p T ( w , k) 



e 2 4\k,t) = -- ^ A A (k) £ A (k) / du; 

4 A=1 J-oc 



ikt 



,—i(u>—k)t 1 



to — k 



t - 



i(uj — fc) 



Akt 



+ 



+ 



where 



2ik i(io — k) 

Pr(w,k) = -Imn r (£j,k) = -^2— it 



+ (OJ — > —u) > + c.c. 



(3.21) 



1 - 



d(k 2 - uo 2 ). 



Potential secular terms (that grow in time) arise at long times from the regions in uj at 
which the denominators in Eq. ( |3.21| ) are resonant. These are extracted in the long time 
limit by using the formulae in an appendix of Ref. [21|. Particularly relevant to the case 
under consideration are the following pl| : 
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r % (1 - cos yt) P (y) 1 =°° J t p(0) + p'(0) [ln(> 1) + 7B ] 
Jo y I 



+ f % \p(v) - p(o) - y p'(o) c(ai - y)] + o (t- 1 ) 



00 dy f smyt\ , . 

— * by 

o y V y / 



o y 

dy 

y 



+ iPvf- [p(y) - p(0) 0(n -y)] + (r l ) , (3.22) 
Jo y v y 



where 7^ = 0.5772157. . . is the Euler-Mascheroni constant and as can be easily shown the 
dependence on /x in the above integrals cancels. For our analysis, the integration variable 
y = oo ± k and p(y) = Pt{^, k)/u;. The terms that grow linearly in time are recognized as 
those emerging in Fermi's golden rule from elementary time-dependent perturbation theory 
in quantum mechanics. We note that p(0) = and p'(0) 7^ 0, thus the first integral above 
gives a contribution to the real part of the mean field (i.e., the amplitude) that features 
a logarithmic secular term, whereas the second integral contributes to the imaginary part 
(i.e., the oscillating phase) with a linear secular term, which as will be recognized below 
determines a perturbative shift of the oscillation frequency. 
Substituting p T (u;,k) into Eq. ( |3.21|) , we obtain 



e 2 4\k,t) = -j2A x (k) 

A=l 



-ikt 



e 2 T 2 
12k 2 



(In 2kt + 7i? - 1) + i 5 k t + O (V 1 ) 



where 



Ren r (a;,k) 



2uo 



e 2 T 2 
12k 



+ c.c. , (3.23) 



(3.24) 



Secular divergences are an ubiquitous feature in the perturbative solution of differential 
equations with oscillatory behaviour. A resummation method that improves the asymptotic 
solution of these differential equations was introduced in Ref. 0] which is interpreted as 
a dynamical renormalization group. This method is very powerful and not only allows a 
consistent resummation of the secular terms in the perturbative expansion but provides 
also a consistent reduction of the dynamics to the slow degrees of freedom as is explicitly 
explained in the work of Chen et. al. and Kunihiro et. al. in |j3~4] |. This method was recently 
generalized to the realm of quantum field theory as a dynamical renormalization group to 
resum the perturbative series for the real time evolution of non-equilibrium expectation 
values 



|2T|,|22|. This generalization of the work of reference [g4[] for differential equations to 
non-equilibrium quantum field theory is a major step since in non-equilibrium systems the 
equations of motion for expectation values are non-local and as shown in detail in this work 
require a resummation of Feynman diagrams. 

While the linear secular terms have a natural interpretation in terms of renormalization 
of the mass (the imaginary part) or a quasiparticle width (the real part) |[21|| , the logarithmic 
secular term found above is akin to those found in Refs. pl| , |22[| that lead to anomalous relax- 
ation. Furthermore the origin of these logarithmic secular terms is similar to the threshold 
infrared divergences and threshold enhancement of the spectral density due to the presence 



of nearby poles 0, 
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Thus following the method presented in Refs. [21,0] we implement the dynamical renor- 
malization group by introducing a (complex) amplitude renormalization factor in the fol- 
lowing manner, 

A x (k) =Z k {r) A x {k,r) , 

where Z k (r) = 1 + e 2 z k (t) + C(e 4 ) is a multiplicative renormalization constant, A\(k, r) 
is the renormalized initial value, and r is an arbitrary renormalization scale at which the 
secular divergences are cancelled [ 2lj^2 |. Choosing 



e 2 T 2 
12k 2 



(In 2kr + ^ E - 1) + i 5 k r , 



we obtain to lowest order in e 2 T 2 /k 2 that the solution of the equation of motion is given by 



a T (k,t) = ]TA(k,T) £ x (k) e 

A=l 

+ nonsecular terms . 



-ikt 



e 2 T 2 



+ c.c. 



(3.25) 



which remains bounded at large times t provided that r is chosen arbitrarily close to t. The 
solution does not depend on the renormalization scale r and this independence leads to the 
dynamical renormalization group equation [ 2l|j22 l, which to this order is given by 



8 e 2 T 2 . c ■ 
d^ + 12k^ + t6k 



A x (k,r) = 0, 



with the solution 



A x (k,r) =A x (k,T ) e 



-iS k (r-T ) 



where ro is the time scale such that this intermediate asymptotic solution is valid and 
physically corresponds to a microscopic scale, i.e, r ~ 1/k. Finally, setting r = t in 
Eq. (|3.25|) we find that ay(k, t) evolves at intermediate asymptotic times t ^> 1/k as 



a r (M) = XM A (k,r ) ^(k) e 



-i(k+5 k )(t-T ) 



A=l 



+ C.C. 



(3.26) 



From Eq. ( |3.24j ) we see that 5 k is consistent with the photon thermal mass m 2 = e 2 T 2 /6 (TT| 
for k 2 3> m 2 . 

As discussed in detail in Ref. [^T| the dynamical renormalization group solution ( |3.26| ) is 
also obtained via the Fourier transform of the renormalization group improved propagator 
in frequency-momentum space, hence the above solution corresponds to a renormalization 
group improved resummation of the self-energy [pl|^2| . 

This novel anomalous power law relaxation of the photon mean field will be confirmed 
below in our study of the kinetics of the photon distribution function in the linearized 
approximation. 

This anomalous power law relaxation is obviously very slow in the semihard regime 
in which the HTL approximation and perturbation theory is valid. At higher orders we 
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expect exponential relaxation due to collisional processes, which emerges from linear secular 
terms 0,|2 in a perturbative solution of the real-time equations of motion. The power 



law relaxation will then compete with the exponential relaxation and we expect a crossover 
time scale at which relaxation will change from a power law to an exponential. Clearly an 
assessment of this time scale requires a detailed calculation of higher order contributions 
which we expect to study in a forthcoming article. 

A related crossover of behavior will be found below for the evolution of the photon 
distribution function and photon production. 



3. Secular terms from the Laplace transform for semihard photons 



We show here how secular terms for semihard photons (eT C fc < T) can be obtained 
from the Laplace transform representation of the photon mean field given by Eqs. (|3.13|) - 
( |3 . 14] ) for large times kt ^> 1. 

The integrand in Eq. ( |3.13| ) can be expanded as follows for semihard momentum, 



1 a < l \ er«fc<r e 2 T 2 



'e 2 T 2 ^ 



1 + 



Inserting this expansion in Eq. ( |3.13| ) gives for a£"*(k, t), 



. cut , 



k,t) 



e 2 T 2 
' 12k 



PV 



+ k dcu 



-iujt 



k UJ^ 



k 2 



a T (k, 0) 



' e 2 r r 
k 2 



1 + 



where PV stands for the principal value prescription. This integral yields for asymptotic 
times to first order in e 2 T 2 /k 2 



(M) 



e 2 T 2 / x e 4 T /r 

cosfct [ln2fct + 7ij] a r (k,0) + C? -, 



12A; 2 



where we used |^T[ 
r°° dy 



(1 - cos yt) p(y) *=°°p(0) [ln(/xt) + 7] + 



y 



V k* 



-\p(y)-p(o) B{j*-v)\ + o0- 
y \t 



The pole contribution (|3.12|) to first order in e 2 T 2 /k 2 reads 

e 2 T 2 



a^ e (k,t) = a T (k,0) 



cos kt 



12 k 2 



t sin kt 



e 4 T 4 ' 
~k± 



+ nonsecular terms + O 



Collecting pole and cut contributions we find 

f e 2 T 2 



& T (k,t) = a T (k, 0) jcos/rf 



12k 2 



[(In 2kt + 'Je) cos kt + t sin kt] 



- e 4 T 4~ 
~kF 



+ nonsecular terms + O 



(3.27) 



This result coincides with the perturbative solution ( |3.20| ), ( [3.231 ) and ( |3.24|) after imposing 
the initial conditions ( |3.5| ) used for the Laplace transform to Eq. ( p.20| ). That is, 

A x (k) = A* x (k) . 

In summary, the Laplace transform solution permits to compute the secular terms as well 
as the perturbative method used in sec. II. A2. 
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B. Quantum kinetics of the photon distribution function 



As mentioned in Sec. II the first step towards a kinetic description is to identify the 
proper degrees of freedom (quasiparticles) and the corresponding microscopic time scale. 
For semihard photons of momenta eT <C k <C T, it is adequate to choose the free photons 
as the quasiparticles with the corresponding microscopic scale ~ 1/k. A kinetic description 
of the non-equilibrium evolution of the distribution function assumes a wide separation 
between the microscopic and the relaxation time scales. In the case under consideration 
the effective small coupling in the semihard momentum limit is e 2 T 2 /k 2 <C 1. Furthermore 
assuming that e C 1 then in this regime both the HTL and the perturbative approximation 
are valid. 

We begin by obtaining the photon number operator from the Heisenberg field operator 
and its conjugate momentum. Write 

A T (x,f) = A r (k,t) e* kx 



P T (x,t) = 



(2tt) 3 / 2 
d 3 k 



Pt(M) e 



ik-x 



with 



Ar(M) 
Pr(M) 



2 

E 

A=l 
2 

*E 

A=l 



2k 

Ob 



a x (k,t) £ x (k) + a[(-k,t) £ x (-k) 



-k,t) £ x (-k)-a x (k,t) £ x (k) 



where a x (k,t) [a x (k,t)] is the annihilation (creation) operator that destroys (creates) a free 
photon of momentum k and polarization A at time t. The polarization-averaged number 
operator iV 7 (k, t) that counts the semihard photons of momentum k is then defined by 



1 

N 7 (k,t) = -J2<(k,t)a x (k,t) 

A A=l 

= l{P T (-k,0-PT(k,0 



A; 2 A T (-k,t) • A T (k,t) 



+ ik[A T (-k, t) ■ P T (k, t) - P T (-k, t) ■ A r (k, *)]} • 

The expectation value of this number operator is interpreted as the number of photons per 
polarization per unit phase space volume 



nj(f) = WM)> = (2tt) 



dA 
d 3 xd 3 k 



where N is the total number of photons per polarization in the plasma. Taking the time 
derivative of iV 7 (k, t) and using the Heisenberg equations of motion, we obtain the following 
expression for the time derivative of the expectation value 



nlit) = lim ^- (-^- - ik 
kV f^tAk \dt' 



d 3 q 



(2tt) 3 / 2 



(^"(-p, t)j ■ A+(k, 0^"(q, 0) + cc , (3.28) 
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where the "+" ( "— " ) superscripts for the fields refer to fields defined in the forward (back- 
ward) time branch in the CTP formulation. We have separated the time arguments and the 
fields on different branches to be able to extract the time derivative from the expectation 
value. 

The nonequilibrium expectation values on the right-hand side (RHS) of Eq. ( |3.28| ) can 
be computed perturbatively in powers of e using the nonequilibrium Feynman rules and 
real-time propagators. 

We assume the initial density matrix is diagonal in the basis of the photon occupation 
numbers with nonequilibrium initial photon populations given by n^(to) and there is no 
initial photon polarization asymmetry. Furthermore we assume that the fermions are in 
thermal equilibrium at a temperature T. 

At 0(e) the RHS of Eq. ( |3.28| ) vanishes identically. This is a consequence of our choice 



of initial density matrix diagonal in the basis of the photon number operator. To 0(e 2 ), we 
obtain 

Kit) = [1 + nl(t )\ r<(t) - nftfo) m) > ( 3 -29) 
where the time-dependent rates are given by 

sin[(a> — k)(t — t )] 



rj(t) = / °° du n* (u, k) 

J —00 



u) — k 



n< (a;, *) = j J [[1 - (k • p) (k • q)]{n F (p) [1 - n F (q)]6(u, -p + q) 

+ [1 - n F (p)]n F (q)5(u +p - q)} + [1 + (k • p)(k • q)] 
x{n F (p)n F (q)5(uj -p - q) + [1 - n F (p)\ [1 - n F (q)]5(uj + p + q)j 



(3.30) 



with p = k + q and TZ^(u,k) is obtained from TZ<(u),k) through the replacement n F <-> 
1 — n F . 

A comment here is in order. As explained above we are focusing on the leading HTL 
approximation, consequently, in obtaining TZ^(u, k) we use the free real-time fermion prop- 
agators which correspond to the hard part of the fermion loop momentum. In general there 
are contributions from the soft fermion loop momentum region which will require to use the 
HTL-resummed fermion propagators for consistency. We will postpone the detailed 
study of the contribution from soft loop momentum to a forthcoming article, focusing here 
on the comparison between the lowest order calculation in real time and the result available 
in the literature for the photon production rate that include HTL corrections. 

Since the fermions are in thermal equilibrium, the Kubo-Martin-Schwinger (KMS) con- 
dition holds: 

n>(u,k)=e^n<(u;,k) , (3.31) 

where (3 = 1/T. It is apparent to see that TZ^ > \uj, k) has a physical interpretation in terms 
of the off-shell photon production (absorption) processes in the plasma. The first two terms 
in TZ^(u,k) describe the process that a fermion (or an anti-fermion) emits a photon, i.e, 
bremsstrahlung from the fermions in the medium, the third term describes annihilation of 
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a fermion pair into a photon, and the fourth terms describes creation of a photon and a 
fermion pair out of the vacuum. The corresponding terms in TZ>(u, k) describe the inverse 
processes. 

As argued, for semihard photons of momenta eT C fc < T, the leading contribution 
of n;{u,k) arises from the hard loop momenta q ~ k. A detailed analysis of TZ^(ui,k) 



along the familiar lines in the hard thermal loop program [pHq,|lT|j shows that in the HTL 
approximation (q ^> k) 



2rp3 / ,2 



^ < (-^)Ihtl = T ^1 v 1-^J 9(e-^). (3.32) 

Thus we recognize that in the HTL limit TZ^^iuj^k) is completely determined to the off- 
shell Landau damping process in which a hard fermion in the plasma emits (absorbs) a 
semihard photon, i.e, bremsstrahlung from the fermions in the medium. 



Dynamical Renormalization Group and the emergence of detailed balance 

We now turn to the kinetics of semihard photons. To obtain a kinetic equation from 
Eq. (|3.29|), we implement the dynamical renormalization group resummation as introduced 



in Refs. pl| , p2| . Direct integration with the initial condition yields 

K(t) = K(t ) + [1 + nl(t )} fdt'T<{t')- nl(t ) fdt'T>{f). (3.33) 

J to J t() 

The integrals that appear in the above expression: 

Vrl(0 = p 1 -<"!(" - W-frll , (3.34) 

t j-oo 1 (lo — ky 

are dominated, in the long time limit, by the regions of uo for which the denominator is 
resonant, i.e, u ~ k. The time dependence in the above integral along with the resonant 
denominator is the familiar form that leads to Fermi's golden rule in elementary time de- 
pendent perturbation theory. In the long time limit t — to ^> 1/k Fermi's golden rule 
approximates the above integrals by 

f df rj(0 « - t )K*(u = k, k) EE . 
Jto 2 

Therefore the rate of photon production, i.e, the coefficient of the linear time dependence 
vanishes at this order because of the vanishing of the imaginary part of the photon self- 
energy on the photon mass shell. However the use of Fermi's golden rule, which is the 
usual approach to extract (time-independent) rates, misses the important off-shell effects 
associated with finite lifetime processes. These can be understood explicitly by using the 
first of formulae (|3~3^ . We find that for t-t ^>l/k 



I if r< (t <) ^ {.„ m - *)] + 7, - 1} + o L^L- 



dt! T> (t') k(t = ]>>1 e 0k f df T< {if) . (3.35) 

to J tg 
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The second line of this equation displays the condition for detailed balance and holds for 
time scales t — t ^> 1/k. We emphasize that this condition is a consequence of the fact that 
the region u ~ k (i.e, the resonant denominator) dominates the long time behavior of the 
integrals in the time-dependent rates Q3.34j ) and the KMS condition Eq. (|3.31|) holds under 
the assumption that the fermions are in thermal equilibrium. 

Thus we obtain one of the important results of this article the time- dependent rates obey 
a detailed balance and that detailed balance emerges in the intermediate asymptotic regime 
t — t ^> 1/k when the secular terms dominate the long time behavior. This is a noteworthy 
result because it clearly states that detailed balance emerges on microscopic time scales 
t — t > 1/k at which the secular terms dominate the integrals but for which a perturbative 
expansion is still valid. Detailed balance then guarantees the existence of an asymptotic 
equilibrium solution which is reached on time scales of the order of the relaxation time. 

The (logarithmic) secular terms in the time-dependent rates ( p.35|) can now be resummed 



using the dynamical renormalization group method by introducing a multiplicative 

renormalization of the distribution function 

n l(to) = 2 k (t ,r) nl(r) , 
Z k (t ,r) = l + e 2 zl}\t ,T)+O(e A ) , 

thus rewriting Eq. ( ^.33[ ) consistently to 0(e 2 ) as 



<{t) = nl(r) + e*4\t ,r) nftr) + [1 + nj(r)] f dt> T<{t>) - tC&t) f dt> T>{t>) + G(e 4 ) . 

Jto Jto 

The renormalization coefficient z^\to, r) is chosen to cancel the secular divergence at a time 
scale t = r. Thus to C(e 2 ), the choice 

e 2 z£\t,r) nl(r) = - [1 + n£(r)] f dt' V<{t') + n£(r) f dt' T>(t') , 

leads to an improved perturbative solution in terms of the "updated" occupation number 



K(r), 



Kit) = nl(r) + [1 + nl{r)\ J* dt' T<(t') - nj(r) [ dt' T>(t>) + 0(e 4 ) , 

which is valid for large times t ^> t provided that r is chosen arbitrarily close to t. A 
change in the time scale r in the integrals is compensated by a change of the K(t) in such 
a manner that n^(t) does not depend on the arbitrary scale r. This independence leads to 
the dynamical renormalization group equation [^l|,0] which consistently to 0(e 2 ) is given 
by 

±-nl{r) = [1 + nl(r)] T<(r) - <(r) T>(r) + G(e 4 ) . (3.36) 



Choosing r to coincide with t in Eq. ( |3.36|) , we obtain the quantum kinetic equation to order 
e 2 given by 

K(t) = [i + K(t)] r<(t) - ni(t) r>(t) . (3.37) 
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For intermediate asymptotic times t — to 3> 1/k at which the logarithmic secular terms 
dominate the integrals for the time-dependent rates and detailed balance emerges, we find 



r£(f) 



fe(t-to)>i e 



fe(*-*o)>i e 



2^3 




6A; 3 t - t 



where the detailed balance relation is explicitly displayed and the terms being neglected are 
oscillatory on time scales ~ 1/k and fall off faster. 

The detailed balance relation between the time- dependent rates ( |3.37| ) guarantees the 
existence of an equilibrium solution of the kinetic equation with the occupation number 

Kit 



oo 



The full solution of this kinetic equation is given by 



7fc (t) = r>(t)-r<(t). 



/ dt'T<(t') e 

Jt 



(3.38) 



However, in order to understand the relaxation to equilibrium of the distribution function, 
we now focus on the linearized relaxation time approximation which describes the approach 
to equilibrium of the distribution function of one k-mode which has been displaced slightly 
off equilibrium while all other modes (and the fermions) are in equilibrium. Thus writing 
n k(^o) = nB{k) + 5ril_(to) with tzb(/c) the Bose-Einstein distribution function, which as shown 
above is an equilibrium solution as a consequence of detailed balance (|3.35|), we obtain 



Shl(t) = -Snl(t) 



For semihard photons of momentum eT <C k <C T, we can simply replace e /3fc — 1 by k/T 
and upon integration we find 



5n1(t) = 5n1(t c 



t-k 

r 



6fc 2 



for k(t - t ) > 1 



(3.39) 



where ro ~ 1/k. 

A noteworthy feature of Eq. ( |3.39|) is that the relaxation of the distribution function 
for semihard photons in the linearized approximation is governed by a power law with an 
anomalous exponent rather than by exponential damping. This result is similar to that 
found in scalar quantum electrodynamics in the Markovian approximation |23| . Comparing 
Eq. (|3.39|) and Eq. ( p. 26 ), we clearly see that the anomalous exponent for the relaxation of 
the photon distribution function in the linearized approximation is twice that for the linear 
relaxation of the mean field ( p.26| ). 

This relationship is well known in the case of exponential relaxation but our analysis 
with the dynamical renormalization group reveals it to be a more robust feature, applying 
just as well to power law relaxation. 

Our approach to derive kinetic equations in field theory is different from the one often 
used in the literature which involves a Wigner transform and assumptions about the separa- 
tion of fast and slow variables P,|35H38|]. Our work differs in many important respects from 
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these other formulations in that it reveals clearly the dynamics of off-shell effects associated 
with non-exponential relaxation. This aspects will acquire phenomenological relevance in 
our study of photon production by off-shell effects in a plasma with a finite lifetime in the 
next section. 



C. Photon production 

An important and phenomenologically relevant byproduct of the study of the kinetics of 
the photon distribution function is to give an assessment of photon production via off-shell 
effects by the lightest (up and down) quarks from the quark-gluon plasma in which quarks 
and gluons (but not photons) are in thermal equilibrium. To the order under consideration 
the only contribution to the photon self-energy is a fermion loop and although we have 
computed it assuming these fermions to be electrons, we can use this result to study photon 
production from thermalized quarks by simply accounting for the electromagnetic charges 
of the up and down quarks and for the number of colors. To this order 0{ot) there are 
no gluon contributions to the photon self-energy and therefore to photon production and 
the rate for photon production obtained from the imaginary part of the photon self-energy 
evaluated on the photon mass shell vanishes. This can be seen directly from Eqs. ( |3.32| ) and 
( |3.34| ) which show explicitly that the usual (time- independent) rate obtained from Fermi's 
golden rule and determined by the imaginary part of the photon self-energy on the photon 
mass shell vanishes. The underlying assumption that motivates the use of Fermi's golden 
rule is that the lifetime of the system is much larger than typical observation times and the 
limit of t — > oo can be taken replacing the oscillatory factors with resonant denominators 
in Eq. ( |3.34|) by a delta function. This obviously extracts the leading time dependence but 
ignores corrections arising from the finite lifetime of the system under consideration. A 
quark-gluon plasma formed in an ultrarelativistic heavy ion collision is estimated to have 
a lifetime ^ 50 fm/c at RHIC energies thus photons are produced by a source of a finite 
lifetime and off-shell effects could lead to considerable photon production during the lifetime 
of the QGP perhaps comparable to those extracted from the usual Fermi's golden rule results 
for infinite lifetime. The focus of this section is precisely the study of this possibility . 

From the expression for the number of photons per phase space volume given by 
Eq. ( |3.38| ) we can obtain the number of photons produced at a given time t per unit phase 
space from an initial (photon) vacuum state, by setting n^(t ) = and including the charge 
of the up and down quarks and the number of colors. Assuming that photons do not ther- 
malize in the quark-gluon plasma, i.e, that their mean free path is much longer than the size 
of the plasma, we neglect the exponentials in the second term of Eq. ( |3.38| ) which are respon- 



sible for building the photon population up to the equilibrium value. To lowest order in the 
electromagnetic fine structure constant we find from Eq. (|3. 38 ) the total number of semihard 



photons at time t per invariant phase space volume summed over the two polarizations given 
by 



, /•* dN tot (t') _ 2k k (t-t )>i 5aT 3 r 



to 



(3.40) 
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This is a noteworthy result and one of the main points of this article, we find photon 
production to lowest (one-loop) order arising solely from off-shell effects: (q — > and 
q — > qj) and annihilation of quarks (qq — > 7). In the HTL limit, valid for semihard photons, 
the leading contribution arises from Landau damping, i.e, from quarks in the medium. 

In the usual approach [|R],[2(J the photon production rate is obtained from the imaginary 
part of the photon self-energy on-shell and receives contributions only at order 0(aa s ) 
and beyond. The lowest order contribution corresponds to a quark loop with a self-energy 
insertion from gluons as well as electromagnetic vertex correction by gluons. As a result the 
imaginary part of the photon self-energy on the mass shell arises from Compton scattering 
(qg — > 57 and qg — > qj) and pair-annihilation (qq — > gj) |19| , |20| . The lowest order 0(aa s ) 
contribution displays an infrared divergence for massless quarks that is screened by the HTL 
resummation leading to the result quoted in Ref. |19| , p0| for the hard photon production 
rate lH 



k 



dN tot (t) 



d 3 x d 3 k dt 



a a 



T 2 e~ k ' T In 



'2.912 k 



+ 1 



on-shell 



18tt 2 a \ g 2 T 

where a s = g 2 /An with g the QCD coupling constant. The total number of photons produced 
per (invariant) phase space at time t is then given by 



I, fdt dNtot{t>) 
t d 3 xd 3 kdt' 



on-shell 



a a s T 2 (t-t )e- k / T In (^^ + 1 



187T 



9 



T 



(3.41) 



which can obviously be interpreted as a linear secular term in the photon distribution func- 
tion, which is the usual situation corresponding to exponential relaxation of the photon 
distribution function P2J as well as for the photon mean field [EHJ] . 



Comparison between on- and off-shell photon production 

The important question to address is how do the on- and off-shell contributions to photon 
production compare? Obviously the answer to this question depends on the time scale 
involved, eventually at long time scales the linear growth of the photon number with time 
from the on-shell contribution will overwhelm the logarithmic growth from the off-shell, but 
at early times the off-shell will dominate. Thus a crossover regime is expected from the off- 
shell to the on-shell dominance to photon production. This crossover time scale will depend 
on the temperature and the energy of the photons as well as on the numerical value of the 
electromagnetic and strong couplings. 

There are two important time scales that are of phenomenological relevance to photon 
production in a quark-gluon plasma. Firstly the assumption that quarks are in thermal 
equilibrium restricts the earliest time scale to be somewhat larger than about 1 — 3 fm/c 
which is the estimate for quark thermalization ||J0]. Secondly the lifetime of the QGP phase 
is estimated to be somewhere in the range 10 — 50 fm/c. Thus the relevant time scale for 
comparison between the off- and on-shell contributions is 1 — 3 fm/ c <t — t < 10 — 50 fm/ c. 



We follow Kapusta et al. [19] in adding 1 to the argument of the logarithm so as to make it agree 
with the numerical result for k ~ T. For details, see discussion below Eq. (41) in Ref fTs|| . 
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Before we begin the comparison, we point out that Eqs. Q3.40| ) and (|3.41|) have different 
regimes of validity in photon momenta. Whereas Eq. ( |3.40| ) is valid in the semihard regime 
eT ^ k ^ T where the HTL approximation is reliable, Eq. ( |3.41| ) is valid in the hard 
regime k > T [|L9,p(J. We have studied numerically the off-shell contribution to photon 



production for hard momentum k ~ T directly in terms of the time-dependent rates given 
by Eqs. (|3.29|) and ( |3.30| ) without using the HTL approximation and compared it to the HTL 
approximation for 1 fm/c < t < 100 fm/c (with to = 0). The result is depicted in Fig. |3|, 
which displays clearly a logarithmic time dependence in both cases for t > 1 —2 fm/c but with 
different slopes, a consequence of the different temperature and momentum dependence of 
the coefficients of the logarithmic time dependence. The HTL approximation overestimates 
the total number of photons by at most a factor two in the relevant regime 1 fm/c < t < 
100 fm/c for hard photons. The reliability of the HTL approximation ( |3.40| ) in the semihard 
momentum region is confirmed by Fig. £|, which compares the number of photons obtained 
from the numerical evaluation with the full time-dependent rates and the result from the 
HTL approximation in the weak coupling and semihard momentum region. 



We can now compare the off-shell contribution to that of Refs. |I9| , |20[| , which is a result of 
on-shell processes and valid in the hard momentum limit. For this we focus on the relevant 
scenario of a quark-gluon plasma with thermalized quarks at a temperature T ~ 200 MeV 
and of lifetime t ~ 10 — 20 fm/c. These are approximately the temperature and temporal 
scales expected to be reached at RHIC. The value of a s at this temperature is not know 
with much certainty but expected to be ^ 1. Following |19[ we choose a s = 0.4 which 
corresponds to g 2 = 5, we find that for hard photons of momenta k ~ T the ratio of photon 
produced by off-shell processes is comparable to that produced by on-shell processes. 

Fig. |5| depicts the ratio of the number of hard (k ~ T) photons produced by off-shell 
processes in the HTL approximation to that from on-shell processes in the regime of lifetimes 
for a quark-gluon plasma phase 5 fm/c < t < 100 fm/c (here we set to — 0). Accounting for 
the overestimate of the HTL approximation for hard photons, we then conclude that off-shell 
processes such as from quarks in the medium (Landau damping) are just as important as 
on-shell processes. This is another of the important results of this work and that cannot 
be obtained from the usual approach to photon production based on the computation of a 
time-independent rate. 

Furthermore, an important point worth emphasizing is that whereas the on-shell calcu- 
lation displays infrared divergences at lowest order, the calculation directly in real time is 
infrared finite because time acts as an infrared cutoff. Thus a calculation of photon produc- 
tion directly in real time does not require an HTL improvement of the fermion propagators 
to cutoff an infrared divergence. However, the region of soft loop momentum does require 
the HTL resummation of one of the fermion propagators and will be studied in a future 
article. 

A thorough study of off-shell effects and their contribution to photon production in a 
wider range of temperatures and momenta including screening corrections to the internal 
quark lines and the study of potential observables will be the subject of a longer study which 
we think is worthwhile on its own and on which we expect to report soon. 
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IV. HARD FERMIONS OUT OF EQUILIBRIUM 



To provide a complete picture of relaxation and non-equilibrium aspects of a hot QED 
plasma in real time, we now focus on a detailed study of relaxation of fermionic mean fields 
(as induced by an adiabatically switched-on Grassmann source) as well as the quantum 
kinetics of the fermion distribution function. In this section we assume that the photons are 
in thermal equilibrium, and since we work to leading order in the HTL approximation we 
can translate the results vis a vis to the case of equilibrated gluons. In particular we seek to 
study the possibility of anomalous relaxation as a result of the emission and absorption of 
magnetic photons. In Refs. it was found that the relaxation of fermionic excitations 

is anomalous and not exponential as a result of the emission and absorption of magnetic 
photons that are only dynamically screened by Landau damping. The study of the real- 
time relaxation of the fermionic mean fields in these references was cast in terms of the 
Bloch-Nordsieck approximation which replaces the gamma matrices by the classical velocity 
of the fermion. In Ref. [|2T| the relaxation of a charged scalar mean field as well as the 
quantum kinetics of the distribution function of charged scalars in scalar electrodynamics 
were studied using the dynamical renormalization group, both the charged scalar mean 
field and the distribution function of charged particles reveal anomalous non-exponential 
relaxation as a consequence of emission and absorption of soft magnetic photons. While 
electric photons (plasmons) are screened by a Debye mass which cuts off their infrared 
contribution, magnetic photons are only dynamically screened by Landau damping and 
their emission and absorption dominates the infrared behavior of the fermion propagator. 

While the dynamical renormalization group has been implemented in scalar theories it 
has not yet been applied to fermionic theories. Thus the purpose of this section is twofold, 

(i) to implement the dynamical renormalization group to study the relaxation and kinetics 
of fermions with a detailed discussion of the technical differences with the bosonic case and 

(ii) to focus on the real-time manifestation of the infrared singularities associated with soft 
magnetic photons. 



A. Relaxation of the fermionic mean field 

The equation of motion for a fermionic mean field is obtained by following the strategy 
described in section II. We begin by writing the fermionic field as 

^(x,t) = V(x,t) + X ± (x,t), 
^(x,t) = (vl/ ± (x,t)). 



Then using the tadpole method JT7|], with an external Grassmann source that is adiabatically 
switched-on from t = — oo and switched-off at t — 0, we find the Dirac equation for the spatial 
Fourier transform of the fermion mean field for t > given by 



/ d 



7- k] V(M) - f dt' £(k,t-f) V(M') =0 , (4.1) 

I J —oo 



where £(k, t — t') is the retarded fermion self-energy. A comment here is in order. To 
facilitate the study and maintain notational simplicity, in obtaining the equation we neglect 
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the contribution from the instantaneous Coulomb interaction which is irrelevant to the 
relaxation of the mean field and only results in a perturbative frequency shift. 

As noted above the relaxation of hard fermions is dominated by the soft photon con- 
tributions, thus in a perturbative expansion one needs to use the HTL-resummed photon 
propagators to account for the screening effects in the medium. To one-loop order but with 
the HTL-resummed photon propagators given in the Appendix, S(k, t — t') reads 



E(k,t-t') 



duo 



i 7 pi(u), k) cos[a;(t — t')} + 7 • kp 2 (uo, k) sin[cj(t — t')} , (4.2) 



where the spectral densities 
d 3 q 



pi(uo,k) = e 



(27t) 3 J-00 



dq [1 + n B (qo) - n F (p)} 



x 



Pt(Qo, q) + « h(Qo, q)\ [5(u - p - q ) + 6(u + p + q )), 
d 3 q 



(4.3) 



/d q r+°° 
77^0 J_ dq [1 + n B (q ) - n F (p)} 



x 



;k-q)(p-q) p T (qo,q) 



k • p 1 

-TT~ pL(lo,q) [6(uj - p - q ) - 5(uj + p + q )} , (4.4) 



with p = k — q. Here pr{qo, q) is the HTL-resummed spectral density for transverse photon 
propagator defined in Eq. (|3.15|) and p~L{qo,q) is the HTL-resummed spectral density for 
longitudinal photon propagator 



Pi{qo, q) = sgn(go) Z L (q) 6[q% - uj 2 L (q)} + (3 L (q , q) 9(q 2 - q%) , 



2rp2 



Pdqo,q) 



go 



6 q 



q^ + ^f 



(2- 2^1n 



q-qo 



+ [ 



7re 2 T 2 901 2 



(4.5) 



where ujl{q) is the plasmon (longitudinal photon) pole and Zi(q) is the corresponding 
residue [PH| . It is worth pointing out that pi( 2 )(u;,k) is an even (odd) function of u, a 
property that will be useful in the following analysis. Furthermore, to establish a connec- 
tion with results in the literature, it proves convenient to introduce the Laplace transform 
of the retarded self-energy E(s,k) just as in Eq. ( |3.6| ) and its analytic continuation E(u;,k) 
as in Eq. ( |3 . 1 1|) which is given by 



E(w,k) 



+00 



dv 



v — to — i0 : 



-7oPi(v, k) + 7-kp 2 (^, k) 



as 



Following the same strategy in the study of the photon mean field, we define a(k,t — t') 





X(k,t-t') = —a(k,t-t') 
and rewrite Eq. (|4.1| ) as an initial value problem 
d 



i l0 - - 7 • k J ^(k, t) + I dt' a(k, t - t') ^(k, t') - a(k, 0) ^(k, t) = 
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with the initial conditions ip(k, 0) = 4>o(k) and ^(k, t < 0) = 0. 

We are now ready to solve the equation of motion by perturbative expansion in powers 
of e 2 just as in the case of the gauge mean field. Let us begin by writing 

^(M) = ^ 0) (k,t) + e 2 ^ (1) (k,t) + 0(e 4 ) , 
a(k, t-t') = e 2 <r«(k, t - + ^(e 4 ), 

we obtain a hierarchy of equations: 
^7o^-7'k^ (0) (k,t) = 0, 



^ - 7 • k j (k, = <j« (k, 0) V (0) (k, f ) - jf (ft' a« (k, t - t') J,® (k, , 



These equations can be solved iteratively by starting from the zeroth-order (free field) solu- 
tion 

V (0) (k,0 = £ [A,(kH(k) e- m + fl.(k)t;.(-k) e ik \ , 

s 

and the retarded Green's function of the unperturbed problem 

0$(M - = {-ilocos[k{t - t')} + 7 • ksin[A;(t - if)]} 0{t - if). 
Here u s {k) and v s (k) are free Dirac spinors that satisfy 

u s (k) ' 



(To - 7 • k) 







*a-k) , 

The solution to the first-order equation is found to be given by 

(k, t) = ^ (1 ' a) (k, t) + ^ (1 ' b) (k, t) , 



where 



e 2 ^ (1 ' a) (k,t) = -tl a (k) t [A(k) « a (k) e~ ikt - B s (k) v s (-k) e 



iktl 



e 2 ^ b \k,t) = - E ———7- 7V,k)<U(k) «,(k) e"* 4 

7T „ J — c 



UJ — k 



t- 



l _ g-i(w— fe)t 

— fc) 



-5 a (k) u 8 (-k) e 



1 _ gi(w-fc)t 

i(a; — fc) 



with 



7"(k) = / 



+°° 1 
— p 2 (^,k) = - Tr 

oo UJ 4 



Re£(0,k)( 7o -7-k) 



7V k) = — [ Pl (c, k) - p 2 (uj, k)] = -— Tr Im£(w, k)( 7o - 7 • k) 
UJ ^uj L 



(4.6) 



(4.7) 



(4.8) 
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In Eq. (|4.6|) secular terms are explicitly linear in time and are purely imaginary, whereas in 
Eq. (4.7) secular terms may arise at long times from the resonant denominators. 

From the form of 7 6 (u;, k) in Eq. ( |4.8| ), we recognize that 7 b (co>, k) evaluated at the fermion 
mass shell u = k is the fermion damping rate computed in perturbation theory ||11|| . It has 
been shown in the literature that due to the emission and absorption of soft quasi-static 
transverse photons which are only dynamically screened by Landau damping, the fermion 
damping rate exhibits infrared divergences near the mass shell in perturbation theory. 

This becomes evident from the following analysis. For soft photons with q Q , q <C T, we 
can replace 

1 + n B (q ) - n F (p) ~ T/q , p~k-qcos9, 
where cos 9 = k • q, thus write 

7V k) = — / 7^ — [(1 - ™s 2 9) M%, Q) + q) 
io J [Ztt) 6 J~q q L 

x5(lu - k + qcosO - q ). (4.9) 

Here we have neglected the subleading pole contributions which corresponding to emission 
and absorption of on-shell photons [|25|| . Recall that for very soft q <C eT the function 
Pt(%, q) /qo is strongly peaked at q = [see Eq. ( |3.17| ) and Fig. [jj, and as q — > it can be 
approximated by 

-Ar(g <g,<z)^^ir as ( 4 - 10 ) 
qo q 2 

The infrared divergences near the fermion mass shell become manifest after substituting 
Eq. ( |4.10|) into Eq. (|4.9|) . The physical origin of the behavior of the function /3t(?o, q)/qo as 
qo ^ q ~* is the absence of a magnetic mass. 

In order to isolate the singular behavior of 7 6 (c<J, k), we follow the steps in Ref. |25[] and 
write 

- MQo,q) = S(q ) [-£ - - 2 j + - Mqo,q) , (4.11) 

qo \q q ^^p) qo 

where ujp = eT/3 is the plasma frequency and ^t(?o>?) denotes the regular part of the 
transverse photon spectral density. Substituting Eq. (|4.11|) into Eq. ( |4.9| ), we can then 
separate 7 fe (o;, k) into an infrared singular part which is logarithmically divergent near the 
fermion mass shell, given by 



, . 7re 2 Tk f d 3 q ( 1 , . 



and a contribution that remains finite and can be expanded near the fermion mass shell, 
given by 

b ne 2 Tk f d 3 q f+i dq 2 
7 {u, k) = / — — / — [/3 L {qo, q) - cos 9 (3 T {qo, q) 

6 iO J [ITX) 6 J-q q 

+ VT(qo,q)] S{u - k + qcos9 - g ) • 
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Using the delta function 8(q cos 9 — q ) to perform the angular integration yields 



7 V, k) u = aT fin + jj + 0(w - fc) 



(4.12) 



where 



J 



a?? / — 

9 go 



So 2 



#&(?o, ?) - -J #r(?o, ?) + fr(?o, g) 



(4.13) 



The above double integral has been computed analytically in Ref. [g5| with the result I 
In 3/2, 



7 V,k)^-aT In 



\uj — k I In 3' 



+ C(cu - fc) 



We are now in position to find the secular terms in -^^(k, t) that emerge in the intermediate 
asymptotic regime. 

The imaginary part of the secular terms in e 2 ^ 1 '^ (k, t) and e 2 ^ 1,b )(k, t) combine into a 
linear secular term given by 

Si, k (t) = T ta6 k t = Tj Tr[ReS(cu,k)( 7o -7-k) 

for the positive (negative) energy spinors respectively, and no further secular terms arise 
from the higher order expansion around the fermion mass shell in Eq. ( |4.12[ ). This purely 
imaginary linear secular term is thus identified with a perturbative shift of the oscillation 
frequency of the mean field [PIJ and is determined by a dispersive integral of the spectral 
densities pi^c^k) which is rather difficult to obtain in closed form but a detailed analysis 
reveals that 8k is finite. 

The real secular terms are more involved. The contribution to 7 b (u;, k) that is finite as 
uj — > k leads to a linear secular term, whereas for the logarithmically divergent contribution 
as uj — > k the following asymptotic result [EHI] becomes useful: 



/•+<*> dy cq ^ ^ _ _ + o , x 

J-oo y 2 



where we have neglected terms that fall off at long times. Thus we find the real part of the 
secular terms to be given by 



SR,k{t) = -aTti \nujpt + j E - 1 + 



hi3' 



Gathering the above results, at large times t 3> 1/up the perturbative solution reads 
V(k,t) = E {t 1 + S h (t)] A s (k) u s (k) e~ ikt + [1 + S* k (t)] B s (k) v s (-k) e ikt } 



+ nonsecular terms 



(4.14) 
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with 



S k (t) 



-a t 



hiupt + 7^ - 1 + 



In 3 s 



T + id k 



Obviously, this perturbative solution breaks down at a time scale r br ~ [aTln(l/a)] 1 . 
To obtain a uniformly valid solution for large times we now implement a resummation of the 



secular terms in the perturbative series via the dynamical renormalization group |2T , 22 



This method is implemented by defining the (complex) amplitude renormalization as 
follows 

A s {k) =Z k {r) A s {k,r) , 
B s (k) 



Z*(t) B s (k,r) 



where Z k (r) = 1 + az k (r) + 0(e A ). The renormalization coefficient az y k '{r) is chosen to 
cancel the secular divergence at a time scale t = r, i.e, «zf (r) = — 5&(r), thus leading to 

V(k,t) = £ {[1 + S k (t) - S k {r)} A s (k,r) u s {k) e~ ikt 

s 

+ [l + S* k (t)-S* k (T)} B s (k,r) v s (-k) e ikt }, 

which remains bounded at large times provided that r is chosen arbitrarily close to t. The 
mean field does not depend on the arbitrary renormalization scale r and this independence 
leads to the dynamical renormalization group equation, which to order 0(a) is given by 



d 


dS k (r) 




dr 


d 


dS* k (r) 




dr 



A(k,r) 



B s (k,r) 



0, 
. 



with solutions 




A s (k,T ) e- taS ^ 
B s (k,r ) e iaSkT 



-aTr[lno;pT+0.12652.. 



where we have replaced 7^ — 1 + ln3/2 = 0.12652 . . ., and r ~ 1/uJp is the time scale such 
that this intermediate asymptotic solution is valid. Hence choosing the renormalization 
point r to coincide with the time t, we find the long time behavior (for ujpt ^> 1) of the 
fermionic mean field for hard momentum is given by 



V(k, t)=J2 [A(k, r ) «.(k) e- iE ^ + B s (k, r ) v s (-k) e iE W 

s 

where E{k) is the pole position shifted by one-loop corrections 



-aTt[\nu) P t+Q. 12652.. 



(4.15) 



E(k) 



k+-Tr 
4 



ReE(w,k)(7o-7-k) 



Eq. ( |4.15[ ) reveals a time scale for the relaxation of the fermionic mean field r IC \ ~ 
[aTln(l/a)]~ 1 , which coincides with the time scale at which the perturbative solution (|4.14| ) 
breaks down . This highlights clearly the nonperturbative nature of relaxation phenomena. 
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This result coincides with that found in Refs. p5| , p6| via the Bloch-Nordsieck approx- 
imation and in scalar quantum electrodynamics 0] using the dynamical renormalization 
group. 

The main purpose of this section was to introduce the dynamical renormalization group 
for fermionic theories. Furthermore, this is another important and relevant example of the 
reliability and consistency of this novel renormalization group applied to real-time nonequi- 
librium phenomena. 



B. Quantum kinetics of the fermion distribution function 

We now study the quantum kinetic equation for the distribution function of hard 
fermions. There has recently been an intense activity to obtain a Boltzmann equation 



for quasiparticles in gauge theories [27|-25] motivated in part by the necessity to obtain a 
consistent description for baryogenesis in non-abelian theories. Boltzmann equations with 
a diagrammatic interpretation were obtained in p7|p9| in which a collision-type kernel de- 
scribes the scattering of hard quasiparticles. In these approaches this collision kernel reveals 
the infrared divergences associated with the emission and absorption of magnetic photons 
(or gluons) and must be cutoff by introducing a relaxation time scale r re i ~ 1/ aT to leading 



logarithmic accuracy [[29 



In a derivation of quantum kinetic equations for charged quasiparticles in scalar quantum 
electrodynamics [^T[] using the dynamical renormalization group, it was understood that the 
origin of these infrared divergences is the implementation of Fermi's golden rule that assumes 
completed collisions and takes the infinite time limit in the collision kernels. The dynamical 
renormalization group leads to quantum kinetic equations in real time in terms of time- 
dependent scattering kernels without any infrared ambiguity. 

In this section we implement this program in spinor QED to derive the quantum kinetic 
equation for hard fermions. There are several important features of our study that must be 
emphasized: (i) as presented in detail in Sec. II, gauge invariance is automatically taken into 
account by working directly with gauge invariant operators, thus the operator that describes 
the number of fermionic quasiparticles is gauge invariant, (ii) a kinetic description relies on 
a separation between the microscopic and the relaxation time scales, this is warranted in a 
strict perturbative regime and applies to hard fermionic quasiparticles, (iii) the dynamical 
renormalization group leads to a quantum kinetic equation in real time without infrared 
divergences since time acts as an infrared cutoff. 

The program begins by expanding the Heisenberg fermion field in terms of creation and 
annihilation operators as 

= e4koc > 

with 

V(k,t) = k(MH(k) + 4(-MK(-k) 

V s 1 
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V> 1 



(M) = */^E [&J(-M)«J(-k) +4(k,tH(k) 



where = Vk 2 + m 2 and 6 s (k, i) [6j(k, t)] is the annihilation (creation) operator that 
destroys (creates) a free fermion of momentum k and spin s at time t. We have retained the 
fermion mass to avoid the subtleties associated with the normalization of massless spinors, 
the massless limit will be taken later. 

The spin-averaged number operator for fermions with momentum k is then given by 

N f (k,t) = Wbt(k,t) b s (k,t) = ^(-M(# + m) 7o V>(k,f) , 

where = (a?k, k). Similarly, the spin-averaged number operator for anti-fermions with 
momentum k is then given by 

N f (k,t) = 5£dJ(k,f) d.(k,f) = ^(k,t)(^-m) 7 o^(-k,t) . 

Taking time derivative of Nf(k,t) and using the Heisenberg equations of motion, we find 

n£(t) = <A>(k,t)> 

= fe£^/ (^{(^" ( - p ' t/)[7 ° A o (-q,*') - 7 • A r (-q,t')] 
x(^ + m) 70 ^ + (k,t))} + c.c. . 



As before here the "+" ( "— " ) superscripts for the fields refer to fields defined in the forward 
(backward) time branch in the CTP formulation. It is straightforward to check that the 
total fermion number (fermions minus antifermions) is conserved. 

In the hard fermion limit ~ T > m we neglect the fermion mass and obtain to 0(e 2 ) 

4W=e 2 / 7^7 [ +OC dq I* dt"{M x {t Q ) cos[{k - p - q )(t" - t )\ 



(2ir) 3 J-oo Jto 
x [p L (q , q) /C+(k, q) + 2 p T (q , q) /C+(k, q)] + Af 2 (t ) 

x cos [(k+p-q )(t" - t ))[p L {q ,q) £r(k,q) + 2p T (q ,q) /C 2 "(k,q)]} , (4.16) 
where M and /C denote respectively the following statistical and kinematic factors 
Ni{t) = [1 - n f k (t)} n f p (t) n B (q ) - n£(t) [1 - n f p (t))[l + n B (q )) , 
Af 2 (t) = [1 - n{(t)][l - n f p (t)} n B (q ) - n£(t) n f p (t) [1 + n B (q )] , 

/Cf(k,q) = l±k-p , /C 2 ± (k,q) = l T k-p± 1 ~[ k ; q ^ . 

1- f(k-q) 

Here Pt{Qo, q) an d Pl{.%, q) are the HTL-resummed spectral densities for the transverse and 
longitudinal photons given by Eqs. Q3.15l ) and (|4.5| ), respectively. 
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In the linearized relaxation time approximation we assume that n£(to) = np{k) + 5n£(to). 
Then upon integrating over t", we obtain 

= -<5«k(«o) fa-R^k) "°l("-*)('-«o)] , (4.17) 

J-oo h{uj — k) 

where 

+ 2p T (<7o,g) ££"(k,q)] <J(w -p - g ) + [pz(?o,g) £f(k,q) 
+ 2p T (g ,g) /CT(k,q)] <5(w + p + g )} . 

Eq. (|4.17|) can be integrated directly to yield 

faiw = ^(ib)(i- r°^R/(«,k) »-«»[(" -*)(«-«■)! I . (4 . 18) 

[ J-00 7r(c<J — Kj^ J 

The time- dependent contribution above is now familiar from the previous discussions, po- 
tential secular terms will emerge at long times from the regions in which the resonant 
denominator vanishes. This is the region near the fermion mass shell u ~ k, where TZf(u, k) 
is dominated by the regions of small q and q , which physically corresponds to emission and 
absorption of soft photons. As before for soft photons with q, go -C T, we can replace 

1 + n B (qo) — n F (p) ~ T/q , p~k — qcos6, JCf ~ 2, K\ ~ 1 — cos 2 6> , 

thus write lZf(u, k) at uj « k as 



d 3 q (i dq 

x5(uj - k + qcosd - g ) • (4.19) 



1Z f (u, k) = 2ne 2 T [(1 - cos 2 0)/? T (?o, g) + (3 L (q , g)] 

J (Zirr J —a Qn 



Note that the double integral in Eq. ( |4.19 ) is exactly the same as that in Eq. (|4.13|) . Thus 
TZf(u},\c) features an infrared divergence near the fermion mass shell as shown in the previous 
subsection. Following the analysis carried out in the preceding subsection we obtain 



K f (u,k) w = k -2aT 



In 



\uj — k\ ln3 



uop 2 



+ 0[(u - kf] . (4.20) 



Substituting Eq. ( |4.20|) into Eq. (|4.18|) , we find the number of fermions at intermediate 
asymptotic times t — to 1/uJp to be given by 

Snl(t) = 5n[{t ){l - 2aT(t - t )[lncop(t - t ) + 0.12652 . . .]} 

+ nonsecular terms. (4-21) 

As in the case of the fermion mean field relaxation [cf. Eq. (|4.14 )], the perturbative solution 
contains a secular term of the form tint. Obviously, the secular term will invalidate the 
perturbative solution at time scales r re i ~ [2aTln(l/o;)] _1 . In the intermediate asymptotic 
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regime 1/k <^t — t <^ r re i, the perturbative expansion can be improved by absorbing the 
contribution of the secular term at a time scale r into a re-definition of the distribution 
function. Hence we apply the dynamical renormalization group method through a renor- 
malization of the distribution function much in the same manner as the renormalization of 
the amplitude in the mean field discussed above, 

<fo£(*o) = Z k (r,t ) 5n f k (T) , 

with 

Z k (r,t ) = l + a z [ k l \rM)+0{a 2 ) . 

The independence of the solution on the time scale r leads to the dynamical renormalization 
group equation which to lowest order in a is given by 

| ^- + 2aT[l + In oj p (t - t ) + 0.12652... ] j 5n[(r) = 0. (4.22) 

which is obviously of the form of a quantum kinetic equation in the linearized relaxation 
time approximation, but with a time- dependent relaxation rate. 

Solving Eq. ( |4.22 ) and choosing r to coincide with t, we obtain the evolution of the 



fermion distribution function at asymptotic times t — to ^> 1 /k to be given by 

Sn[(t) = Snl(t ) exp{-2aT(t - t )[lnu P (t - t ) + 0.12652 ...]}. (4.23) 



Comparing Eq. ( |4.23|) with Eq. ( 4.15 ) we find that the anomalous exponent that describes 



the relaxation of the fermion distribution function in the linear approximation is twice that 
for the linear relaxation of the mean field. A similar relation is obtained between the damping 
rate for the single quasiparticle relaxation and the relaxation rate of the distribution function 
in the case of time- independent rates and true exponential relaxation |i"Tfl . The dynamical 



renormalization group reveals this to be a generic feature even with time-dependent rates. 



V. CONCLUSIONS AND FURTHER QUESTIONS 

The goals of this article are the study of nonequilibrium effects in hot QED plasmas 
directly in real time. The focus is a systematic study of relaxation of mean fields as well 
as the distribution functions for photons and fermions. In particular the application of the 
dynamical renormalization group method to study anomalous relaxation as a consequence of 
the exchange of soft photons. To begin with, we have cast our study solely in terms of gauge 
invariant quantities, this can be done in the abelian theory in a straightforward manner, 
avoiding potential ambiguities associated with gauge invariance. The relaxation of photon 
mean fields revealed important features: for soft momentum k ^ eT mean fields that are 
prepared by adiabatically switching-on an external source there is exponential relaxation 
towards the oscillatory behavior dominated by the transverse photon pole. The source that 
induces the mean field in this case has a Fourier transform that is singular at zero frequency 
and excites the resonance in the Landau damping region near zero frequency for the soft 
photon mean field. Sources that have a regular Fourier transform would not lead to the 
exponential relaxation. For semihard momentum eT <C k <C T in principle both the HTL 
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and the perturbative approximations are valid, however the spectral density for photons 
becomes sharply peaked at the edge of the Landau damping continuum consistent with the 
fact that the photon pole becomes perturbatively close to the Landau damping cut. This 
enhancement of the spectral density near the bare photon mass shell results in the breakdown 
of perturbation theory at large times kt ^> 1. The dynamical renormalization group provides 
a consistent resummation of the photon self-energy in real time and leads to anomalous 
relaxation of the mean field, given by Eq. ( |3.26| ). Clearly, higher order terms beyond the 
HTL approximation will include collisional contributions leading perhaps to exponential 
relaxation. We then expect a crossover between the anomalous power law obtained in 
lowest order and the exponential relaxation from collisional processes, the crossover time 
scale will depend on the details of the different contributions and requires a study beyond 
that presented here. 

The dynamical renormalization group provides a consistent and systematic framework 
to obtain quantum kinetic equations directly in real time from the microscopic field the- 
ory This method allows to extract information that is not available in the usual 
kinetic description in terms of time-independent collision kernels obtained under the as- 
sumption of completed collisions which only include on-shell processes. The dynamical 
renormalization group approach to quantum kinetics consistently includes off-shell processes 
and accounts for time- dependent collisional kernels. This is important and potentially phe- 
nomenologically relevant in the case of the quark-gluon plasma which has a finite lifetime. 

We have implemented this approach to obtain the quantum kinetic equation for the 
distribution function of semihard photons in the case of thermalized fermions in lowest order 
in the HTL approximation. This equation features time-dependent rates and we established 
that detailed balance, a consequence of fermions being in thermal equilibrium, emerges 
on microscopic time scales. The linearization of the kinetic equation describes relaxation 
towards equilibrium with an anomalous exponent, twice as large as that of the photon mean 
field in the semihard case. 

The kinetic equation for semihard photons allows us to study photon production from 
a thermalized quark-gluon plasma by simply replacing the fermions by two flavors of light 
quarks. To leading order in the HTL approximation photon production arises from off- 
shell bremsstrahlung (q — > and q — > qj), and the total number of photons produced 
at time t is given by Eq. ( |3.40| ). We study the hard region k ~ T numerically and find 
that for T ~ 200 MeV the number of photons produced by these off-shell processes form 
a thermalized QGP of lifetime 10 fm/c < t < 50 fm/c is comparable to that by on-shell 
processes estimated in Refs. [|T9| , pO[| [See Fig. |!|. 

The relaxation of the fermion mean field for hard momentum is studied with the dynam- 
ical renormalization group, we find an anomalous exponential relaxation which confirms the 
results of Refs. [pB1 , |26| where the Bloch-Nordsieck approximation was used. 

We then obtain the quantum kinetic equation for the distribution function of hard 
fermions assuming that photons are in thermal equilibrium. The collisional kernel is time- 
dependent and infrared finite. The linearized kinetic equation describes approach to equi- 
librium with an anomalous exponential relaxation, which is twice that of the fermion mean 
field. 

An important payoff of this approach to quantum kinetics is that it bypasses the assump- 
tion of completed collisions which lead to collisional kernels obtained by Fermi's golden rule 
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and only describe on-shell processes as in the usual kinetic approach, which in the case under 
consideration leads to infrared divergent collisional kernels [pSR . 

Further questions. Perhaps the most phenomenologically pressing aspect that requires 
further and deeper study is the photon production by off-shell processes in a thermalized 
quark-gluon plasma. This is important in view of the fact that RHIC will begin physics runs 
very soon hence an assessment of potential experimental electromagnetic signatures is very 
relevant. The next step in the study of photon production must be to include hard thermal 
loop correction to the quark propagators, and to study the spectral distribution of the 
photons produced as a function of momentum, temperature and lifetime of the thermalized 
QGP. Furthermore, our calculations just as those in Refs. [19|]2(| do not account for the 
expansion of the plasma, thus an important next step is to couple the quantum kinetic 
equation to the hydrodynamic equations of evolution of the QGP, this is required to obtain 
the photon distribution function integrated over the space-time volume of the QGP. Work 
on these aspects is currently underway. 

Another important extension of this work is to obtain the relaxation of non-abelian gauge 
mean fields beyond the HTL approximation both in the soft and semihard regime. Similarly 
for the fermionic mean fields as well as the kinetics of the fermion distribution function. Of 
particular importance in this regard would be to include the nonequilibrium evolution of the 
distribution function of the gauge fields in the collision kernels for the fermionic distribution 
function. Obviously this will require a consistent analysis of the relevant time scales. These 
aspects bear on the physics of baryogenesis and we expect to report on some of our studies 
on this issues soon. 
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APPENDIX: REAL-TIME PROPAGATORS 

In this appendix we summarize the various CTP real-time propagators used in this 
article. 

(i) The free fermion propagators (with zero chemical potential) are defined by 

(vr(x,t)*V,0> = ij -0^sz\t,t')^-*'\ 

S+ + (t, t') = S>(t, t')9(t - + S<(t, t ! )6{t' - t), 
S k -(t,f) = S>(t,t')9(t' -t) + S<(t,t')6{t - 0, 
S^(t,t') = St(t,t'), (Al) 
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where a, b = ± and the fermion Wightman functions read 



S>(t,t" 



2u k 



S<(t,t" 



+ lo ( i K- m ) lo n F (k) e^"')}, 
-{(?C + m) n F (k) e -^(*-*') 



2cj, 



+ 7o(^-m) 7 o[l-^(A;)]e^('-'')}, 



(A2) 



with K = (a>k, k), cuk = Vk 2 + m 2 and np(k) the Fermi-Dirac distribution, 
(ii) The free (transverse) photon propagators are defined by 



(A^(K,t)A^',t')) 



d 3 q 
(2^)3 



^(^0^'(q)e iq - 



iq-(x— x') 



= $T, q (t, t'W -t) + G^ g (t, t')9(t - 0, 



(A3) 



where (q) = 5 y — is the transverse projector and the photon Wightman func- 
tions read 



2? 



-iq(t-t') 



+ [l + n B (g)]e' 



iq(t-t') 



(A4) 



with nsiq) the Bose-Einstein distribution. 



(iii) The HTL-resummed transverse photon propagators are obtained from the free ones 
3) by using the following Wightman functions 



ST, g (t, = i / rf <?o Pr(<7o, (?) [1 + n s (g )] e-^ ^'), 



(A5) 



where is the resumed spectral density, 
(iv) The HTL-resummed longitudinal photon propagators are given by [2 

ffiftf ) = 4*(* - + GUt,t')9(t - + GZAtJW - 1), 



(A6) 
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with the Wightman functions expressed in terms of the resummed spectral density pL 
as 

GtJtJ) =-%J dq p L (q ,q) [1 + n B (q )} e'^^, 
GL, q W) =~ij dqop L (qo,q)n B (q )e- i ^ t - t '\ (A7) 

The free real-time propagators used in studying quantum kinetics are obtained from 
lT|)-( |A4| ) with the equilibrium Fermi-Dirac and Bose-Einstein distributions replacing by 
the initial nonequilibrium ones. Finally, we note that the HTL-resummed photon propa- 
gators are only valid for photons in thermal equilibrium because in deriving their spectral 
representations (KB) and (|A7| ) used has be made of the KMS condition [21,22 . 
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FIG. 1. (3t(u, k)/uj vs lo for ultrasoft photon of momentum k/eT = 0.1. 
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FIG. 2. (3t{w, k)/uo vs u for semihard photon of momentum eT/k = 0.1. 
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FIG. 3. Off-shell contribution to nH£) vs t for hard photons k ~ T, T = 200 MeV and 
a = 1/137. The solid line is the numerical evaluation using the full time dependent rates, the 
dashed line is the HTL approximation. 
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FIG. 4. The number of semihard photons is plotted as a function of time for e = 0.01, 
T = 200 MeV and eT/k = 0.1. The solid line is the numerical evaluation using the full time 
dependent rates, the dashed line is the HTL approximation. 



47 



10 




1 1 1 1 1 1 1 1 1 1 1 

20 40 60 80 100 

t (fm/c) 

FIG. 5. The ratio of the number of hard photons produced by off-shell process (q — > q^f and 
q — > (77) to that produced by on-shell processes (qg — ► 97 and — > 57) as a function of QGP 
lifetime t for a = 1/137, a s = 0.4 and T = 200 MeV. 
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